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BIOENFORMATICALLY DETECTABLE GROUP OF NOVEL 
VIRAL REGULATORY GENES AND USES THEREOF 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 
The present invention relates to a grovq> of bioinfonnatically detectable novel viral 
genes, here identified as Viral Genomic Address Messenger or VGAM genes, which are believed 
to he related to flie micro RNA (miRNA) groiq) of genes. 



DESCRiraON OF PRIOR ART 

MIR genes are regulatory genes encoding MicroRNA's (miRNA), short ~22nt 
non-coding RNA's, found in a wide range of species, believed to function as specific gene 
translation repressors, sometimes involved in cell-differentiation. Some 110 human MIR genes 
have been detected by laboratory means. Over the past 6 months, the need for computerized 
detection of MIR genes has been recognized, and several informatic detection engines have been 
reported (Lim, 2 003; Grad, 2 003; Lai, 2 003). Collectively these informatic d election engines 
found 38 more human MIR genes which were later confirmed in zebrafish, 14 human MIR genes 
which were confirmed in human, and 55 postulated human MIR genes which could not be 
confirmed by laboratory (Lim, 2003). Extensive efforts to identify novel MIR genes using 
conventional biological detection techniques such as massive cloning and sequencing efforts, and 
several bioinfonnatic detection attempts, have led leading researchers in tiie field to the 
conclusion that the total number of human MIR genes is between 200 to 255 (Lau, 2003; Lim 
2003 Science; Lim, 2003 Genes Dev). Recent studies postulate that the number of MIR genes 
may be higher (Grad, 2003; Krichevsky, 2003). 

The ability to detect novel MIR genes is limited by the methodologies used to 
detect such genes. All MIR genes identified so far either present a visibly discemable whole 
body phenotype, as do Lin-4 and Let-7 (Wightman ,1993; Reinhart , 2000), or produce suflScient 
quantities of RNA so as to be detected by the standard molecular biological techniques. 

Initial studies reportmg ME. genes (Bartel, 2001; Tuschl, 2001) discovered 93 
MIR genes in several species, by sequencing a limited number of clones (300 by Bartel and 100 
by Tuschl) of small segments (i.e. size fractionated) RNA- MiRNA encoded by MIR genes 
1 
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detected in these studies therefore, represent the more prevalent among the miRNA gene family, 
and can not be much rarer than 1% of all small '-20nt-long RNA segments. 

Cuirent methodology has therefore been unable to detect micro SNA genes (MIR 
genes) which either do not present a visually discemable whole body phenotype, or are rare (e.g. 
rarer than 0.1% of all size firactiomated ~20nt-long RNA segments expressed in the tissues 
examined), and therefore do not produce significant enough quaittities of RNA so as to be 
detected by standard biological techniques. To date, miKNA have not been detected in viruses. 



BRIEF DESCRIPTION OF SEQUENCE USTING, LARGE TABLES AND COMPUEBR 

PROGRAM LISTING 

Sequence listing, large tables related to sequence listing, and computer program 
listing are filed under section 8 01 (a)(i) on an electronic medium in computer readable form, 
attached to title present invention, and are hereby incorporated by reference. Said sequence 
listing, large tables related to sequence listing, and computer program listing are submitted on a 
CD-ROM (Operating system: MS-Windows), entitled SEQUENCE LISTING AND LARGE 
TABLES, containing files the names and sizes of which are as follows: 

Sequence listing comprising 424,595 geDomic sequences, is filed under section 
801 (a)(i) on an electromc medium in computer readable form, attached to the present invention, 
and is hereby mcorporated by reference. Said sequence listing is contained in a self extracting 
compressed file named SBQ_IIST.EXE (8,389KB). Compressed file contains I file named 
SEQ_LIST.TXT (63,778KB). 

Large tables relating to genomic s equences are stored in 7 self extracting files, 
each comprising a respective one of the following table files: TABLE1.TXT (229KB); 
TABLE2.TXT (878KB); TABLBS.TXT (147KB); TABLE4.TXT (353,556KB); TABLE5.TXT 
(868,334KB); TABLE6.TXT (208,163KB); and TABLE7.TXT (5ZKB). 

It is appreciated that the nucleotide is represented as "T' in the sequences 
incorporated in the enclosed large tables. 

Computer program Usting of a computer program constructed and operative in 
accordance with a preferred embodiment of the present invention is enclosed on an electronic 
medium in computer readable form, and is hereby incorporated by reference. The computer 
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program, listing is contaioed in a self extracting compressed file named COMPUTER 
PROGRAM LISTING.EXE (100KB). Compressed file contains 7 files, the name and sizes of 
which are as foUows: AUXILARY_FILES.TXT (117K3); BIlSIDING_SrrE_SCORING.TXT 
(17KB); EDrr_DISTANCE.TXT (104KB); FIRST-K.TXT (48K); HAIRPIN_PREDICnO.TXT 
(4TKB); TWOJPHASED_PREDICT0R.TXT (74KB); and 

TWO_PHASED_SIDE_SELECTOR.TXr (4KB). 



3 
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SUMMARY OF TEB INVENTION 
The present invention relates to a novel group of bioinfotmatically detectable, 
viral regulatory SNA genes, which repress expression of host target host genes, by means of 
complementary hybridization to binding sites in untranslated regions of these host target host 
genes. It is believed that this novel group of viral genes represent a pervasive viral mechanism of 
attacking hosts, and that therefore knowledge of this novel group of viral genes may be useful in 
preventing and treadng viral diseases. 

In various preferred embodiments, the present invention seeks to provide improved 
method and system for detection and prevention of vnal disease, which is mediated by this group 
of novel viral genes. 

Accordinglly, the invention provides several substantially pure nucleic acids (e.g., 
genomic nucleic acid, cDNA or synthetic nucleic acid) each encoding a novel viral gene of the 
VGAM group of gene, vectors comprising the nucleic acids , probes comprising the nucleic acids 
, a method and system for selectively modulating translation of known "target" genes utilizing the 
vectors, anld a method and system for detecting expression of known "target" genes utilizing the 
probe. 

By "substantially pure nucleic add" is meant nucleic acid that is &ee of the genes 
which, in the naturally-occurring genome ofthe organism from which thenucleic acidofthe 
invention is derived, flank the genes discovered and isolated by the present invention. The term 
&erefore includes, for example, a recombinant nucleic acid which is incorporated into a vector, 
into an autonomously replicating plasmid or virus, or into the genomic nucleic acid of a 
prokaryote or eukaiyote at a site other than its natual site; or which exists as a separate molecule 
(e.g., a cDNA or a genomic or cDNA fiagnent produced by PGR or restriction endonuclease 
digestion) independent of other sequences. It also includes a recombinant nucleic acid which is 
part of a hybrid gene encoding additional polypeptide sequence. 

"Inhibiting translation" is defined as the ability to prevent syutfaesis of a specific 
protem encoded by a respective gene, by means of inhibiting the translation of the mKNA of this 
gene. "Translation inhibiter site" is defined as the minimal nucleic acid sequence sufQcient to 
inhibit translation. 

There is thus provided in accordance with a prefeired embodiment of the present 
invention a bioinformatically detectable novel viral gene encoding substantially pure nucleic acid 
wherdn: SNA encoded by the bioinformatically detectable novel viral gene is about 18 to about 
4 
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24 nucleotides in length, and originates ftom an RNA precursor, which RNA precursor is about 
50 to about 120 nucleotides in length, a nucleotide sequence of a first half of the RNA precursor 
is a partial inversed-reversed sequence of a nucleotide sequence of a second half thereof a 
nucleotide sequence of the KNA encoded by the novel viral gene is a partial inversed-reversed 
sequence of a nucleotide sequence of a binding site associated with at least one host target gene, 
and a fimction of the novel viral gesne is bioinformatically deducible. 

There is further provided in accordance with another preferred embodiment of the 
present invention a method for anti-viral treatment comprising neutralizing said KNA. 

Further in accordance with a preferred embodiment of the present invention the 
neutralizing comprises: synthesizing a complementary nucleic acid molecule, a nucleic sequence 
of which complementary nucleic acid molecule is a partial inversed-reversed sequence of said 
RNA, and transfecting host cells with the complementary nucleic acid molecule, thereby 
oomplementaiily binding said RNA. 

Further in accordance with a preferred embodiment of the present invention tiie 
neutralizing comprises immunologically neutralizing. 

There is still further provided in accordance with another preferred embodiment of the 
present invention a bioinformatically detectable novel viral gene encoding substantially pure 
nucleic acid wherein: RNA encoded by the bioinfonnatically detectable novel viral gene includes 
a plurality of RNA sections, each of the RNA sections being about 50 to about 120 nucleotides in 
length, and including an RNA segment, which RNA segment is about 18 to about 24 nucleotides 
in length, a nucleotide sequence of a first half of each of the RNA sections encoded by the novel 
viral gene is a partial inversed-reversed sequence of nucleotide sequence of a second half thereoj^ 
a nucleotide sequence of each of the RNA segments encoded by the novel viral gene is a partial 
inversed-reversed sequence of the nucleotide sequence of a binding site associated with at least 
one target host gene, and a fimction of the novel viral gene is bioinformatically deducible fi:om 
the following data elements: the nucleotide sequence of the RNA encoded by the novel viral 
gene, a nucleotide sequence of the at least one target host gene^ and fimction of the at least one 
target host gene. 

Further in accordance with a preferred embodiment of the present invention the 
fimction of the novel viral gene is bioinibtmatically deducible fiom the following data elements: 
the nucleotide sequence of the RNA encoded by the bioinformatically detectable novel viral 
gene, a nucleotide sequence of the at least one target host gene, and a fimction of the at least one 
target host gene. 

5 
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Still furfher in accordance with a preferred embodiment of the present invention the 
RNA encoded by the novel viral gene complementarily binds the binding site associated with the 
at least one target host gene, thereby modulatiiig expression of the at least one target host gene. 

Additianally in accordance with a preferred embodiment of the present invention the 
binding site associated with at least one target host gene is located in an imtraoslated region of 
RNA encoded by the at least one target host gene. 

Moreover in accordance with a preferred embodiment of the preseat invention the 
function of the novel viral gene is selective inhibition of translation of the at least one target host 
gene, which selective inhibition inchides complementary hybridization of the BNA encoded by 
the novel viral gene to the binding site. 

Further m accordance with a preferred embodiment of the present invention the 
invention includes a vector inchiding the DNA. 

Still further in accordance with a preferred embodiment of the present invention the 
invention includes a method of selectively inhibiting translation of at least one gene, including 
introducing the vector. 

Moreover in accordance with a preferred embodiment of the present invention the 
introducing includes utilizing RNAi pathway. 

Additionally in accordance with a preferred embodiment of the present invention the 
invention includes a gene expression inhibition system including: the vector, and a vector 
inserter, functional to insert the vector into a cell, thereby selectively inhibiting translation of at 
least one gene. 

Further in accordance with a preferred embodiment of the present invention the 
invention includes a probe including the DNA. ° 

Still further in accordance with a preferred embodiment of the present irrvention the 
invention includes a method of selectively detecting expression of at least one gene, including 
using the probe. 

Additionally in accordance with a preferred enibodimeirt of the present invention the 
invention includes a gene expression detection system includiing: the probe, and a gene 
expression detector functional to selectively detect expression of at least one gene. 

Further in accordance with a preferred embodiment of the present invention the 
invention includes an anti-viral substance capable of neutralizing the RNA. 

Still further in accordance with a preferred embodiment of the present invention the 
neutralizing includes complementarily binding the KNA. 

6 
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Additionally in accordance with a preferred embodiment of the present invention the 
neutralizing includes immunologically neutralizing. 

Moreover in accordance with a preferred embodiment of the present invention the 
invention includes a method for anti-viral treatment including neutralizing the KNA. 

Further in accordance with a preferred embodiment of the present invention the 
neutralizing includes: S3mthesizing a complementary nucleic add molecule, a nucleic sequence of 
which complementary nucleic acid molecule is a partial inversed-reversed sequence of the KNA, 
and transfecting host cells with the complementary nucleic acid molecule, thereby 
complementarily binding the KNA. 

Still further in accordance with a preferred embodiment of the present invention the 
neutralizing includes immunologically neutraliziiig. 

BRIEF DESCRIPTION OF DRAWINGS 
Fig. 1 is a simplified diagram illustrating a mode by which viral genes of a novel 
groiQ) of viral genes of the present invention, modulate expression of known host target genes; 

Hg. 2 is a simplified block diagram illustrating a bioinformatic gene detection 
system capable of detecting genes of the novel groupofviral genesof the present invention, 
which system is constructed and operative in accordance with a preferred embodiment of the 
present invention; 

Fig. 3 is a simplified flowchart illustrating operation of a mechanism for training 
of a computer system to recognize the novel viral genes of the present invention, which 
mechanism is constructed and operative in accordance with a preferred embodiment of the 
present invention; 

Fig. 4A is a simplified block diagram of a non-coding genomic sequence detector 
constructed and operative in accordance with a preferred embodiment of the present invention; 

Fig. 4B is a sunpUfied flowchart illustrating operation of a non-coding genomic 
sequence detector constructed and operative in accordance with a preferred embodiment of the 
presMit inveation; 

Fig. 5A is a sin^lified block diagram of a hairpin detector constructed and 
operative in accordance with a preferred embodiment of the present invention; 
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Fig. 5B is a simplified flowchart illustraling operation of a haiipin detector 
constructed and operative in accordance with a preferred embodiment of the present invention; 

Fig. 6A is a simplified block diagram of a dica:-cut location detector constructed 
and operative in accordance with a preferred embodiment of flie present invention; 

Fig. 6B is a simplified flowchart illustrating training of a dicer- cut location 
detector constructed and operative in accordance with a preferred embodiment of the present 
invention; 

Fig. 6C is a simplified fiowchait illustrating operation of a dicer-cut location 
detector constructed and operative in accordance with a prefisired embodiment of the present 
invention; 

Fig. 7A is a simplified block diagram of a target-gene binding-site detector 
constructed and operative in accordance with a preferred embodiment of the present invention; 

Fig. 7B is a simpliSed flowchart illustrating operation of a target-gene bindingnsite 
detector constructed and operative in accordance with a preferred embodimeait of the present 
invention; 

Fig. 8 is a simplified flowchart illustrating operation of a function & utility 
analyzer constructed and operative in accordance with a prefened embodiment of the present 
invention; 

Fig. 9 is a simplified diagram describing a novel bioinfi>imatically detected group 
of viral regulatory genes, referred to here as Viral Genomic Record (VGR) genes, each of which 
encodes an 'operon-like' cluster of novel iniRNA-Iike viral genes, which in turn modulate 
expression of one or more host target g^es; 

Fig. 10 is a block diagram illustrating different utilities of novel viral genes and 
novel operon-like viral genes, both of the present invention; 

Figs. 1 1 A and 1 IB are simplified diagrams, which when taken together illustrate a 
mode of gene therapy applicable to novel viral genes of the present invention; 

Fig. 12 is a table summarizing laboratory validation results which validate eflScacy 
of a bioinformatic gene detection system constructed and operative in accordance with a 
preferred embodiment of the present invention; 
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Fig. 13 is a picture of laboratoiy results validating fixe ejqiression of 37 novel 
human genes detected by a bioinfonnatic gene detection eaguiG constructed and operative in 
accordance with a prefened embodiment of the present invention, thereby validating the efficacy 
of the gene detection engine of the present inventioi:^ 

Fig. 14A is a schematic representation of an 'operon like' cluster of novel human 
gene hairpin sequences detected bioinformatically by a bioinformatic gene detection engine 
constructed and operative in accordance with a preferred embodiment of the present invention, 
and non-GAM hairpin useful as negative controls thereto; 

Fig. 14B is a schematic lepresentatian of secondary folding of hairpins of the 
operon-like cluster of Fig. 14A4 

Fig. 14C is a picture of laboratory results demonstrating expression of novel 
human genes of Figs. 14A and 14B, and lack of expression of the neg&tive controls, thereby 
validating efficacy of bioinformatic detection of GAM genes and GR genes of the present 
iavention, by a bioinformatic gene detection engine constructed and operative in accordance with 
a preferred embodiment of the present invention; 

Fig. 15A is an annotated sequence of EST72223 comprising known human 
miRNA gene MIR98 and novel human gene GAM25, both detected by the gene detection system 
of the present invention; and 

Figs. ISB, 15C and 15D are pictures of laboratory results demonstrating 
laboratory confirmation of ejqnression of known human gene MUt^S and of novel 
bioinformaticaUy detected human gene GAM25 respectively, both of Fig. 15A, tbus validating 
the bioinformatic gene detection system of the present invention. 

Fig, 16 presents pictures of laboratory results demonstrating laboratory 
confirmation of 'didng' of four novel bioinformatically detected HIVl VGAMs into their 
correspondmg mature genes, herein designated VGAM2032.2 (Fig. 16B), VGAM3249.1 
(Fig.l6C), VGAM507.2 (Fig.l6D) and VGAM1016.2 (Fig.l6E). 

Fig. 17 presents pictures of laboratory results demonstrating laboratory 
confirmation of expression of two novel bioinformatically detected Vaccinia VGAM precursors, 
herein designated VGAM224, and VGAM3184. 
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DETAILED DESCRIPTION OF DRAWINGS 

Reference is now made to Fig. 1, which is a simplified diagram describing each of 
a plurality of novel 'bioinfonnatically detected viral genes of the present invention, refeired to 
here as Viral Genomic Address Messenger (VGAM) genes, which modulates expression of 
respective target genes thereof, the function and utility of which host target genes is known in the 
art. 

VGAM is a novel bioinfonnatically detected regulatory, non protein coding, micro 
RJSfA (miRNA) viral gene. The method by which VGAM is detected is described hereinabove 
with reference to Figs. 1-8. 

VGAM GENE is gene contained in the virus genome and TARGET GENE is a 
human gene contained in the human genome. 

VGAM GENE encodes a VGAM PRECURSOR RNA. Similar to other miKNA 
genes, and unlike most ordinary genes, VGAM PRECURSOR RNA does not encode a protein. 

VGAM PRECURSOR RNA jfolds onto itself, forming VGAM FOLDED 
PRECURSOR RNA, which has a two-dimensional ^hairpin structure'. As is well known in the 
art, this 'hairpin structure', is typical of RNA encoded by miRNA genes, and is due to the feet 
that the nucleotide sequence of the first half of the RNA raicoded by a miRNA gene is an 
accurately or partially inversed reversed sequence of the nucleotide sequence of the second half 
thereof. By inversed reversed is meant a sequence which is reversed and wherein each nucleotide 
is replaced by a complementary nucleotide, as is well known in the art (e.g ATGGC is the 
reverse complementary sequence of GCCAT). 

An enzyme complex designated DICER COMPLEX, 'dices' the VGAM 
FOLDED PRECURSOR RNA into VGAM RNA, a single stranded -22 nt long RNA segment. 
As is known in the art, 'dicing' of a hairpin structured RNA precursor product into a short ~22nt 
RNA segment is catalyzed by an enzyme complex comprising an enzyme called Dicer together 
with other necessary proteins. 

TARGET GENE encodes a corresponding messaiger RNA, VGAM TARGET 
RNA. VGAM TARGET RNA comprises three regions, as is typical of mRNA of a protein 
coding gene: a 5' untranslated region, a protein coding region and a 3' untranslated region, 
designated 5'UTR, PROTEIN CODING and 3'UTR respectively. 
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VGAM RNA binds complementarily to one or more target binding sites located in. 
untranslated regions of VGAM TARGET RNA. This complementary blading is due to the fact 
that the micleotide sequence of VGAM RNA is a partial or accurate inversed reversed sequence 
ofthenucleotide sequence of each ofthe host target binding sites. As an illustration. Fig. 1 
shows three such target binding sites, designated BINDING SITE I, BINDING SITE H and 
BINDING SITE HI respectively. It is appreciated that the number of host target binding sites 
sihown in Fig. 1 is meant as an illustration only, and is not meant to be limiting VGAM RNA 
may have a di^arent number of host target binding sites in untranslated regions of a VGAM 
TARGET RNA It is further appreciated that while Fig. 1 depicts host target binding sites in the 
3'"UTR region, this is meant as an example only, these target binding sites may be located in the 
3'UTR region, the 5'UTR region, or in both 3'UTR and 5'UTR regions. 

The c omplementary b inding o f VGAM RNA tot arget b inding s ites o n VGAM 
TARGET RNA, such as BINDING SITE I, BINDING SITE n and BINDING SITE m, inhibits 
translation of VGAM TARGET RNA into VGAM TARGET PROTEIN. VGAM TARGET 
PROTEIN is therefore outlined by a broken line. 

It is appreciated that TARGET GENE in jEact represents a plurality of VGAM 
target genes. The mRNA of each one of this plurality of VGAM target genes comprises one or 
more target binding sites, each having a nucleotide sequence which is at least partly 
con^ilementary to VGAM RNA, and which when bound by VGAM RNA causes inhibition of 
translation of respective one or more VGAM host target proteins. 

It is ftirther appreciated by one skilled in the art that the mode of transIationaJ 
inhibition illustrated by Fig. 1 with specific reference to translational inhibition exerted by 
VGAM GENE on one or more TARGET GENE, is in fact common to other known miRNA 
genes, as is well known in the art. 

Nucleotide sequences of each of a plurality of VGAM GENEs described by Fig. 1 
and their respective genomic source and chromosomal location are further described hereinbelow 
with reference to Table 1, hereby incorporated by reference. 

Nucleotide sequences of VGAM PRECURSOR RNA, and a schematic 
representation of a predicted secondary folding of VGAM FOLDED PRECURSOR RNA, of 
each of a pluraUty of VGAM GENEs described by Fig. 1 are flirlher described hereinbelow with 
reference to Table 2, hereby incorporated by reference. 
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Nucleotide sequences of a 'diced' VGAM RNA of each, of a plurality of VGAM 
GENEs described by Fig. 1 aire further described hereinbelow with reference to Table 3, hereby 
incorporated by reference. 

Nucleotide sequences of host target binding sites, such as BINDING SITB-I, 
BINDING SUE-n and BINDING STTE-III of Fig. 1, found on VGAM TARGET RNA, of each 
of a plurality of VGAM GENEs described by Fig. 1, and schematic representation of the 
complementarity of each of these host target binding sites to each of a plurality of VGAM RNA 
described by Fig. 1 are described hereinbelow with reference to Table 4, hereby incorporated by 
reference. 

It is appreciated that specific functions and accordingly utUities of each of a 
plurality of VGAM GENEs described by Fig. 1 correlate with, and may be deduced from, the 
identity of the TARGET GENEs that each of said plurality of VGAM GENEs binds and inhibits, 
and the function of each of said TARGET GENEs, as elaborated hereinbelow with reference to 
Table 5, hereby incorporated by reference. 

Studies establishing known functions of each of a plurality of TARGET GENEs of 
VGAM GENEs of Fig. 1, and correlation of said each of a plurality of TARGET GENEs to 
known diseases are listed in Table 6, and are hereby incorporated by reference. 

The present invention discloses a novel group of genes, the VGAM genes, 
belonging to the miRNA genes group, and for which a specific complementary binding has been 
determined. 

Reference is now made to Fig. 2 which is a simplified block diagram illustrating a 
bioinformatic gene detection system CE^able of detecting genes of the novel group of genes of the 
present invention, which system is constructed and operative in accordance with a preferred 
embodiment of the present invention. 

A c enterpiece o f the p resent i nvention i s a b ioinformatic gene d etection engine 
100, which is a prefBrred implementation of a mechanism capable of bioinformaticaUy detecting 
genes of the novel group of genes of the present invention. 

The fiinction of the bioinformatic gene detection engine 100 is as follows: it 
receives three types of input, expressed RNA data 102, sequenced DNA data 104, and protein 
fimction data 106, performs a complex process of analysis of this data as elaborated below, and 
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based on this analysis produces ou^ut of a bioinfoimatically detected group of novel genes 
designated 108. 

Expressed KMA data 102 comprises published expressed sequence tags (EST) 
data, published mRNA data, as well as other sources of published RNA data. Sequenced DNA 
data 104 comprises alphanumeiic data describing sequenced genomic data, which preferably 
includes annotation data such as location of known protein coding regions relative to the 
sequenced data. Protein function data 106 comprises scientific pubUcations reporting studies 
which elucidated physiological function known proteins, and their connection, involvement and 
possible utiUty in treattneat and diagnosis o f v arious diseases. E xpressed RNA data 1 02 and 
sequenced DNA data 1 04 may preferably be obtained from data published by the National Center 
for Bioinfonnatics (NCBI) at the National histitute of Health (NH) (Jenuth, 2000), as well as 
from various other published data sources. Protein function data 106 may preferably be obtained 
from any one of numerous relevant published data sources, such as the Onhne Mendelian 
Inherited DiseaseInMan< OMtMCTM)) d atabase d eveloped b y Jo hn H opkins U niversity, and 
also pubUshed by NCBI (2000). 

Prior to actual detection of bioinfonnatically detected novel genes 108 by the 
bioinformatic gene detection engine 100, a process of bioinformatic gene detection aigine 
training & validation designated 110 takes place. This process uses the known miRNA genes as 
a training set (some 200 such genes have been found to date using biological laboratory means), 
to train the bioinformatic gene detection engine 100 to bioinformatically recognize miRNA-like 
genes, and their respective potential host target binding sites. Bioinformatic gene detection 
engine training & validation 1 10 is further described hereinbelow with reference to Fig. 3. 

The bioinformatic gene detection engine 100 comprises several modules which are 
preferably activated sequentially, and are described as follows: 

A non-coding genomic sequence detector 112 operative to bioinformatically detect 
non-protein coding genomic sequences. The non-coding genomic sequence detector 112 is 
further described herein below with reference to Figs. 4A and 4B. 

A hairpin detector 114 operative to bioinformatically detect genomic 'haiipin- 
shaped' sequences, similar to VGAM FOLDED PRECURSOR of Fig. 1. The hairpui detector 
1 14 is further described herein below with reference to Figs. 5A and 5B. 
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A dicer-cut location detector 116 operative to bioinformatically detect the location 
on a haiipin shaped sequence which is enzymatically cut by DICER COMPLEX of Fig. 1. The 
dicer-cut location detector 1 16 is further described herein below with reference to Fig. 6A. 

A target-gene binding-site detector 118 operative to bioinformatically detect host 
target genes having binding sites, the nucleotide sequence of which is partially complementary to 
that of a given genomic sequence, such as a sequence cut by DICER COMPLEX of Fig. 1. The 
target-gene binding-site detector 118 is fiirther described hereinbelow with reference to Figs. 7A 
andTB. 

A function & utility analyzer 120 operative to analyze function and utility of 
target genes, in order to identify host target genes which have a signijScant clinical function and 
utility. The function & utility analyzer 120 is fiirfher described hereinbelow with reference to 
Fig. 8. 

ELardware implementation of the bioinfbrmatic gene detection engine 100 is 
important, since significant computing power is preferably required in order to perform the 
computation of bioinformatic gene detection engine 100 in reasonable time and cost. For 
example, it is estimated that a using a powerful 8-processor server (e.g. DELL POWERBDGE 
(TM) 8450, 8 XEON (TM) 550MHz processors, 8 GB RAM), over 6 years (!) of computing time 
are required to detect all MIR genes in the human EST data, together with their respective 
binding sites. Various computer hardware and software configurations may be utihzed in order 
to address this computation challenge, as is known in the art. A preferred embodiment of the 
present mvention may preferably comprise a hardware configuration, comprismg a cluster of one 
hundred PCs (PENTIUM (TM) IV, 1.7GHz, with 40GB storage each), connected by Ethemet to 
12 servers (2-CPU, XEON (TM) 1.2-2.2GHz, with ~200GB storage each), combmed with an 8- 
processor server (8-CPU, Xeon SSOMhz w/ 8GB RAM) coimected via 2 HBA fiber-channels to 
an EMC CLARIION (TM) lOO-disks, 3.6 Terabyte storage device. A prefisTred embodiment of 
the present invention may also prefiarably comprise a software configuration which utihzes a 
commercial database software program, such as MICROSOFT (TM) SQL Server 2000. Using 
such preferred hardware and software configuration, may reduce computing time required to 
detect all MIR genes in the human EST data, and their respective binding sites, fiom 6 years to 
45 days. It is ^preciated that the above mentioned hardware configuration is not meant to be 
limitmg, and is given as an illustration only. The present invention may be implemented in a 
wide variety of hardware and software configurations. 
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The present invention discloses 1560 novel genes of the VGAM group of genes, 
which have been detected bioinfonnatically, as described hereinbelow with reference to Table 1 
througli Table 6, and 205 novel genes of the GR group of genes, which have hem detected 
bioinfbimaticaJly, as described hereinbelow with reference to Table 7. Laboratoiy confiimation 
of 37 bioinformatically predicted genes of the human GAM and GR group of genes, and several 
bioinfonnatically predicted genes of the VGAM and VGR group of genes, is described 
herembelow with reference to Figs. 13 through 17. 

Reference is now made to Fig. 3 which is a simplified flowchart illustrating 
operation of a mechanism for training a computer system to recognize the novel genes of the 
present invendon. This mechanism is a preferred implementation of the bioinformatic gene 
detection engine training & validation 1 10 described hereinabove with reference to Fig. 9. 

BIOINFORMATIC GENE DETECTION ENGINE TRAINING & 
VALIDATION 110 of Fig. 2 begins by training the bioinfonnatic gene detection engine to 
recognize Icnown miRNA genes, as d esignated by numeral 1 22. T his training step comprises 
HAIRPIN DETECTOR TRAESIING & VALIDATION 124, further described hereinbelow with 
reference to Fig. 5A, DICER-CUT LOCATION DETECTOR TRAINING & VALIDATION 
126, further described hereinbelow with reference to Fig. 6A and 6B, and TARGET-GENE 
BINDING-SrrE DETECTOR TEIAINING & VALIDATION 128, fiirther described hereinbelow 
with reference to Fig. 7A. 

Next, the BIOINFORMATIC GENE DETECTION ENGINE 100 is used to 
bioinfonnatically detect sample novel genes, as designated by numeral 130. Examples of sample 
novel genes thus detected are described hereinbelow with reference to Fig. 12. 

Finally, wet lab experiments are preferably conducted in order to validate 
expression and preferably function of the sample novel genes detected by the BIOINFORMATIC 
GENE DETECTION ENGINE 100 in the previous step. An example of wet-lab validation of 
the above mentioned sample novel gene bioinfonnatically detected by the system is described 
hereinbelow with reference to Figs. 13 throu^ 17. 

Reference is now made to Fig. 4A which is a simplified block diagram of a 
preferred unpleanentatian of the NON-CODING GENOMIC SEQUENCE DETECTOR 112 
described hereinabove with reference to Fig. 2. The NON-PROTEIN CODING GENOMIC 
SEQUENCE DETECTOR 112 of Fig. 2 preferably receives as input at least two types of 



15 



wo 2004/048511 PCT/IL2003/000998 
pubUshed genomic data: EXPRESSED RNA DATA 102 and SEQUENCED DNA DATA 104. 
The EXPRESSED RNA DATA can include, among others, EST data, EST clusters data, EST 
g^ome alignment data and mRNA data. Sources for EXPRESSED RNA DATA 102 include 
NCBI dbEST, NCBI UniGene clusters and mapping data, and TIGR goie indices. SEQUENCED 
DNA DATA 104 includes both sequence data (FASTA format files), and features annotation 
(GenBank file format) mainly fiom NCBI database. After its initial training, indicated by numeral 
134, and based on the above mentioned ii^ut data, the NON-PROTEIN CODING GENOMIC 
SEQUENCE DETECTOR 112 produces as output a plurality of NON-PROTEIN CODING 
GENOMIC SEQUENCES 136. Preferred operation of the NON-PROTEIN CODING 
GENOMIC SEQUENCE DETECTOR 112 is described hereinbelow with reference to Fig. 4B. 

Reference is now made to Fig. 4B which is a simplified flowchart illustrating a 
preferred operation of the NON-CODING GENOMIC SEQUENCE DETECTOR 112 of Fig. 2. 
Detection of NGN PROTEIN CODING GENOMIC SEQUENCES 136, generally preferably 
progresses in one of the following two paths: 

A first path for detecting NON-PROTEIN CODING GENOMIC SEQUENCES 
136 begins by receiving a plurality of Imown RNA sequences, such as EST data. Each RNA 
sequence is first compared to all known protein-coding sequences, in order to select only those 
RNA sequences which are non-protein coding, i.e. intergenic or mtronic. This can preferably be 
perfbimed by sequence comparison of the RNA sequence to known protein coding sequences, 
using one of many alignment algorithms known in the art, such as BLAST. This sequence 
comparison to the DNA preferably also provides the localization of the KNfA sequence on the 
DNA. 

Alternatively, selection of non-protein coding RNA sequences and their 
localization to the DNA can be performed by using publicly available EST clusters data and 
genomic mapping databases, such as UNIGENE database published by NCBI or TIGR database, 
in order map expressed RNA sequences to DNA sequences encoding them, to find the right 
orientation of EST sequences, and to exclude ESTs which map to protein coding DNA regions, 
as is well known in the art. Public databases, such as TIGR, may also be used to map an EST to a 
cluster of ESTs, assumed to be expressed as one piece, and is known in the art as Tentative 
Human Consensus. Publicly available genome annotation databases, such as NCBFs 
GENBANK, may also be used to deduce expressed intronic sequences. 
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Optionally, an attempt may be made to 'expand' the non-protein KNA sequences 
thus found, by searching for transcription start and eud signals, upstream and doiwnstreain of 
location of the RNA on the DNA. respectively, as is well known in the art 

A second path for detecting non-protein coding geaaomic sequences starts by 
receiving DNA sequences. The DNA sequences are parsed into non protein coding sequences, 
based on published DNA annotation data, by extracting those DNA sequences which are between 
known protein coding sequences. Next, transcription start and end signals are sought. If such 
signals are found, and depending on their 'strength', probable expressed non-protein coding 
genomic sequences are yielded. Such approach is especially useful for identifying novel GAM 
genes which are found in proximity to other known miRNA genes, or other wet-lab validated 
GAM genes. Since, as described hereinbelow with reference to Fig. 9, VGAM genes are 
frequently found in clusters, therefore sequences near a known nuRNA are more likely to contain 
novel genes. Optionally, sequence orthology, i.e. sequences conservation in an evolutionary 
related species, may be used to select genomic sequences having higher probability of containing 
expressed novel VGAM genes. 

It is appreciated that in the present invention the bicdnfonnatics gene detection engine 100 
utilize the input genomic sequences, without filtering protein coding regions detected by the non- 
coding genomic sequence detector 112. 

Reference is now made to Fig. 5A which is a simplified block diagram of a 
preferred implemoitation of the HAIRPIN DETECTOR 114 described hereinabove with 
reference to Fig. 2. 

The goal of the HAIRPIN DETECTOR 1 14 is to detect 'hairpin' shaped genomic 
sequences, similar to those of known miRNA genes. As mentioned hereinabove with reference to 
Fig. 1, a 'hairpin' genomic sequence refers to a genomic sequence which 'folds onto itself 
forming a hairpin like shape, due to the fact that nucleotide sequence of the first half of the 
nucleotide sequence is an accurate or partial complementary sequence of the nucleotide sequence 
of its second half. 

The HAIRPIN DETECTOR 114 of Fig. 2 receives as input a plurality of NON- 
PROTEIN CODING GENOMIC SEQUENCES 136 of Fig. 4A. After a phase of HAIRPIN 
DETECTOR TRAINING & VALIDATION 1 24 o f Fig. 3 , the HAIRPIN DETECTOR 1 14 i s 
opearative to detect and ou^ut 'hairpin shaped' sequences, which are found in the input NON- 
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PROTEHSr CODING GENOMIC SEQUENCES 138. The haiipiu shaped sequences detected by 
the HAIRPIN DETECTOR 114 are designated HAIEIPINS ON GENOMIC SEQUENCES 138. 
Prefened operation of the HAIRPIN DETECTOR 1 14 is descaibed hereinbelow with reference to 
Fig. 5B. 

The phase of HAIRPIN DETECTOR TRAINING & VALIDATION 124 is an 
iterative process of applying the HAIRPIN DETECTOR 114 to known hairpin shaped miRNA 
genes, cahT^rating the HAIRPIN DETECTOR 1 14 such that it identifies the training set of known 
hairpins, as weU as sequences which are similar thereto. In a preferred embodiment of the 
present invention, THE HAIRPIN DETECTOR TRAINING & VALIDATION 124 trains and 
validates each of the steps of operation of the HAIRPIN DETECTOR 114, which steps are 
described hereinbelow with refsrence to Fig. 5B. 

The hairpin detector training and validation 124 prefisrably uses two sets of data: a 
training set of known miRNA genes, such as 440 m iRNA genes of H. sapiens, M. musculus, C. 
elegans, C. Brigssae and D. Melanogaster, annotated in RFAM database (Griffiths- Jones 2003), 
and a large 'background set' of haiipins found in expressed non-protetn coding genomic 
sequences, sucii as a set of 21,985 hairpins found in Tentative Human Concensus (THC) 
sequences in TIGR database. The 'background set' is expected to comprise some vahd, 
previously undetected miRNA hairpins, and many hairpins which are not miRNA hairpins. 

In order to validate the performance of the HAIRPIN DETECTOR 114, a 
validation method is preferably used, which validation method is a variation on the k-fold cross 
validation method ^tchell, 1997). This preferred validation, method is devised to better cope 
with the nature of the training set, which includes large famiUes of similar and even identical 
miRNAs, The training set is preferably first divided into clusters of miRNAs such that any two 
miRKAs that belong to difEbrent clusters have an Edit Distance score (see Algorithms and 
Strings, Dan Gusfield, Cambridge University Press, 1997) of at least D=3, i.e. they differ by at 
least 3 editing operations. Next, the group of clusters is preferably divided into k sets. Then 
standard k-fold cross validation is preferably peiformed on this gmvip of clusters, preferably 
using 1^, such that the members of each cluster are all in the training set or in ibs test set. It is 
appreciated that without the prior clustering, standard cross validation methods results in much 
higiher performance of the predictors due to the redundancy of training examples, within the 
genome of a species and across genomes of different species. 
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Iq a preferred embodiment of the present invention, using the abovementioned 
validation method, the efBcacy of tiie HAIRPIN DETECTOR 114 is indeed validated: for 
example, when a similarity threshold is chosen such that 90% of the published miRNA-precuisor 
hairpins are successfully predicted, only 7.6% of the 21,985 background hairpins are predicted to 
be miRNA-precursors, some of which may indeed be previously unknown miSNA. precursors. 

Bfifarence is now made to Fig. 5B which is a simplified flowchart illustrating a 
preferred operation of the HAIRPIN DETECTOR 1 14 of Fig. 2. 

A hairpin structure is a secondary structure, resulting from the nucleotide 
sequence pattern: the nucleotide sequence of the first half of the hairpin is a partial or accurate 
inversed reversed sequence of the nucleotide sequence of the second half thereof. Various 
methodologies are known in the art for prediction of secondary and tertiary hairpin structures, 
based on given nucleotide sequences. 

In a preferred embodiment of the present invention, the HAIRPIN DETECTOR 
114 initially calculates possible secondary structure folding patterns of a given one of the non- 
protein coding genomic sequences 136 and the respective energy of each of these possible 
secondary folding patterns, preferably using a secondary structure folding algorithm based on 
free-energy minimization, sudi as the MFOLD algorithm (Mathews et al., 1999), as is well 
known in the art. 

Next, the HAIRPIN DETECTOR 114 analyzes the results of the secondary 
structure folding, in order to determine the presence, and location of hairpin folding structures. A 
secondary structure folding algorithm, such as MFOLD algorithm, typically provides as output a 
listing of the base-pairing of the folded shape, i.e. a listing of each pair of connected nucleotides 
in the sequence. The goal of this second step is to asses this base-pairing listing, m order to 
determine if it describes a hairpin type bonding pattern. Pre&rably, each of the sequences that is 
determined to describe a hairpin structure is folded separately in order to determine its exact 
folding pattern and free-energy. 

The HAIRPIN DETECTOR 114 then assess those haiipin structures found by the 
previous step, comparing them to hairpins of known miRNA genes, using various characteristic 
hairpin features such as length of the haiipin and of its loop, free-energy and thermodynamic 
stability, amount and type of mismatched nucleotides, existence of sequence repeat-elements. 
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Only hairpins tiiat bear statistically significant resemblance to the training set of known miRNA 
hairpins, according to the abovementioned parameters are accepted. 

In a prefened embodiment of the present mvention, similarity to the training set of 
known miRNA hairpins is determined using a 'similarity score' which is calculated using a 
weighted sum of terms, where each term is a jftinction of one of the abovementioned hairpin 
features, and the parameters of each function are learned from the set of known hairpins, as 
described hereinabove with reference to hairpin detector training & validation 124. The weight 
of each temi in the similarity score is optimized so as to achieve maximal separation between the 
distribution of similarity scores of hairpins which have been validated as miRNA-precursor 
hairpins, and the distribution of similarity scores of hairpins detected in the 'background set' 
mentioned hereinabove with reference to Fig. SB, many of which are expected not to be mdRNA- 
precursor hairpins. 

In another pire&ired embodiment of the present inventiQn, the abovemeaitioned 
DETERMINE IF SIMILAR TO KNOWN HAIRPIN-GENES step may may preferably be spht 
mto two stages. The first stage is a pennissive filter &at implements a simplified scoring method, 
based on a subset of the haupin features described hereinabove, such as minimal length and 
maximal free energy. The second stage is more stringent, and a JEull calculation of the weighted 
sum of terms described heremabove is performed. This second stage may preferably be 
performed only on the subset of hairpins that survived prior filtering stages of the hairpin- 
detector 114. 

Lastly, the HAIRPIN DETECTOR 114 attempts to select those hairpin structures 
which are as theimodynamically stable as the hairpins of known miRNA genes. This may 
preferably be achieved in various manners. A preferred embodiment of the present invention 
utilizes the following methodology preferably comprising three logical steps: 

First, the HAIRPIN DETECTOR 114 attempts to group potential hairpins into 
•families' of closely related hairpins. As is known in the art, a free-energy calculation algorithm, 
typically provides multiple 'versions' each describing a different possible secondary structure 
folding pattern for the given genomic sequence, and tiie free energy of such possible folding, 
The HAIRPIN DETECTOR 1 14 therefore preferably assesses all hairpuis found in each of the 
'versions', grouping hairpins which appear in different versions, but which share near identical 
locations ioto a common 'family' of hairpins. For example, all hairpins in diJBEerent versions, the 



20 



wo 2004/048511 PCT/IL2003/000998 
center of which hairpiiis is within 7 nucleotides of each other may preferably be grouped to a 
single '£amily'. Ilaiipias may also he grouped to a single '&tmly' if the sequences of one or more 
hairpins are identical to, or are subsequences of, the sequence of another hairpin. 

Next, haiipin 'families' are assessed, in order to select only those families which 
represent hairpins that are as stable as those of known miKNA hairpins. Preferably only famihes 
which are represented in a majority of the secondary structure folding versions, such as at least in 
65% or 80% or 100% of the secondary structure folding versions, are considered stable. 

Finally, an attempt is made to select the most suitable haiipin from each selected 
farmly. For example, a haiipin which appears in more versions than other hairpins, and in 
versions the free-energy of which is lower, may be preferred. 

In another preferred embodiment of the present invention, hairpins with homology 
to other species, and clusters of thermodynamicaUy stable hairpin are fiirther favored. 

Reference is now made to Fig. 6A which is a simplified block diagram of a 
preferred implementation of the DICER-CUT LOCATION DETECTOR 116 described 
hereinabove with reference to Fig. 2. 

The goal of the DICER-CUT LOCATION DETECTOR 116 is to detect the 
location in which DICER COMPLEX of Fig. 1, comprising the enzyme Dicer, would 'dice' the 
given haiipin sequence, similar to VGAM FOLDED PRECURSOR KNA, yielding VGAM RNA 
both of Fig. 1. 

The DICER-CUT LOCATION DETECTOR 116 of Fig. 2 therefore receives as 
input a plurality of HAIRPINS ON GENOMIC SEQUENCES 138 of Fig. 5 A, which were 
calculated by the previous step, and after a phase of DICER-CUT LOCATION DETECTOR 
TRAINING & VALIDATION 126, is operative to detect a respective plurality of DICER-CUT 
SEQUENCES FROM HAIRPINS 140, one for each hairpin. 

in a preferred embodiment of the present invention, tiie DICER-CUT LOCATION 
DETECTOR 116 preferably uses standard machine learning techniques such as K nearest- 
neighbors, Bayesian networks and Support Vector Machines (SVM), trained on known dicer-cut 
locations of known miRNA genes in order to predict dicer-cut locations of novel VGAM genes. 
The DICER-CUT LOCATION DETECTOR TRAINING & VALIDATION 126 is fiarther 
described hereinbelow with reference to Fig. 6B. 
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Reference is now made to Fig. 6B which is a shiiphfied flowchart illustratfaig a 
prefeared implementation of DICBR-CUT LCX^ATION DETECTOR TRAINING & 
VALIDATION 126 of Fig. 3. 

The general goal of the DICER-CUT LOCATION DETECTOR TRAINING & 
VALIDATION 126 is to analyze known haiipin shaped miRNA-precuisors and their respective 
dicer-cut miRNA, in order to determine a common pattern to the dicer-cut location of known 
miRNA genes. Once such a common pattern is deduced, it may preferably be used by the 
DICER-CUT LOCATION DETECTOR 116, in detecting the predicted DICER-CUT 
SEQUENCES FROM HAIRPINS 140, fiom the respective HAIEIPINS ON GENOMIC 
SEQUENCES 138, all of Fig. 6A. 

First, the dicer-cut location of all known miRNA genes is ohtamed and studied,' so 
as to train the DICER-CUT LOCATION DETECTOR 116: for each of the known miRNA, the 
location of the miRNA relative to its hairpin-shaped miRNA-precursor is noted. 

The 5' and 3' ends of the dicer-cut location of each of the known roiRNA genes is 
represented relative to the respective miRNA precursor hairpin, as well as to the nucleotides in 
each location along the hairpin. Frequency and identity of nucleotides and of nucleotide-paiiing, 
and position of nucleotides and nucleotide pairing relative to the dicer-cut location in the known 
miRNA precursor hairpins is analyzed and modeled. In a preferred embodiment of the present 
mvention, features learned j&ora pubUshed miRNAs include: distance from hairpin's loop, 
nucleotide content, positional distribution of nucleotides and mismatched-nucleotides, and 
synunetry of mismatched-nucleotides. 

Different techniques are well known in the art of machine learning for analysis of 
existing pattern from a given 'training set' of examples, which techniques are then capable, to a 
certain degree, to detect similar pattems in other, previously unseen examples. Such machine 
learning techniques include, but are not limited to neural networks, Bayesian networks. Support 
Vector Machines (SVM), Genetic Algorithms, Markovian modeling. Maximum Liklyhood 
modehng. Nearest Neighbor algorithms, Decision trees and other techniques, as is well knovm in 
the art. 

The DICER-CUT LOCATION DETECTOR 116 preferably uses such standard 
machine learning techniques to predict either the 5' end or both the 5' and 3' ends of the miRNA 
excised, or *diced' by the Dicer enzyme from the miRNA haiipin shaped precursor, based on 
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known pairs of miRNA-precursors and their respective resisting miRNAs. The nucleotide 
sequences of 440 published miRNA and their corresponding hairpin p recursors are preferably 
used for training and evaluation of the dicer-cut location detector module. 

Using the abovementioned training set, machine learning predictors, such as a 
Support Vector Machine (SVM) predictor, are implemented, which predictors test every possible 
nucleotide on a hairpin as a candidate for being the 5' end or the 3' end of a miRNA. Other 
machine learning predictors include predictors based on Nearest Neighbor, Bayesian modeling, 
and K-nearest-neigihbor algorithms. The training set of the published miRNA precursor 
sequences is preferably used for training multiple separate classifiers or predictors, each of which 
produces a model for the 5' or 3' end of a miRNA relative to its hairpin precursor. The models 
take into account various miRNA properties such as the distance of the respective (3' or 5') end 
of the miRNA from the hairpin's loop, the nucleotides at its vicinity and the local 'bulge' (i.e. 
base-pair mismatch) structure. 

Petfoimance of the resulting predictors, evaluated on the abovementioned 
validation set of 440 published miRNAs using k-fold cross validation (Mitchell, 1997) with k = 
3, is found to be as follows: in 70% of known miRNAs S'-eod location is correctly determined by 
an SVM predictor within up to 2 nucleotides; a Nearest Neighbor (EDIT DISTANCE) predictor 
achieves 53% accuracy (233/440); a Two-Phased predictor that uses Baysian modeling (TWO 
PHASED) achieves 79% accuracy (348/440), when only the first phase is used, and 63% 
(277/440) when both phases are used; a K-nearest-neighbor predictor (FIRST-K) achieves 61% 
accuracy (268/440). The accuracies of all predictors are considerably higher on top scoring 
subsets of pubUshed miRNA. 

Finally, in order to validate the efficacy and accuracy of the dicer-cut location 
detector 116, a sample of novel genes detected thereby is preferably selected, and validated by 
wet lab. Laboratory results validating the efficacy of the dicer-cut location detector 116 are 
described hereinbebw with reference to Figs. 12 through 15D. 

Reference is now made to Fig. 6C which is a simplified flowchart iUustrating 
operation of DICER-CUT LOCATION DETECTOR 116 of Fig. 2, constructed and operative in 
accordance with a preferred embodiment of the present invention. 

The DICER CUT LOCATION DETECTOR 116 is a machine learning computer 
program module, which is trained on recognizing dicer-cut location of known miRNA genes, and 
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based on this training, is operable to detect dicer cut location of novel VGAM FOLDED 
PRECURSOR RNA. In a preferred embodiment of the present invention, the dicer-cut location 
module preferably utilizes machine learning algorithms, such as Support Vector Machine (SVM), 
Bayesiaa modeling. Nearest Neighbors, and K-nearest-neighbor, as is well known in the art. 

When assessing a novel VGAM precursor, all 19-24 nucleotide long se^ents 
comprised in the VGAM precursor are initiaUy considered as 'potential VGAMs', since the 
dicer-cut location is initially unknown. 

For each such potential VGAM, its 5' end, or its 5' and 3' ends are scored by two 
or more recognition classifiers or predictors. 

In a preferred embodiment of the present invention, the DICER-CUT LOCATION 
DETECTOR 116 preferably uses a Support Vector Machine predictor trained on features such as 
distance from haiipiu's loop, nucleotide content, positional distribution of nucleotides and 
mismatched-nucleotides, and symmetry of mismatched-nucleotides. 

In another preferred embodiment of the present invention, the DICER-CUT 
LOCATION DETECTOR 116 preferably uses an 'EDIT DISTANCE' predictor, which seeks 
sequences that are similar to those of pubhshed miRNAs, utilizing the Nearest Neighbor 
algorithm, where the similarity metric between two sequences is a variant of the edit distance 
algorithm (Algorithms and Strings, Dan Gusflsld, Cambridge University Press, 1997). This 
predictor is based on Ihe observation that miRNAs tend to form clusters (Dostie, 2003), the 
members of which show mark^ sequence similarity to each other. 

In yet another preferred embodiment of the present invention, the DICER-CUT 
LOCATION DETECTOR 116 preferably uses a 'TWO PHASED' predictor, which predicts the 
dicer-cut location in two distinct phases: (a) selecting the double-stranded segment of the hairpin 
comprising the miKNfA by naive Bayesian modeling (Mitchell, 1997), and (b) detecting which 
strand contains the miRNA by either naive or by K-nearest-neighbor modeling. The latter is a 
variant o f t he 'FIRST-K' p redictor d escribed h erein b elow, w ifh p arameters o ptimized f or t his 
specific task. The TWO PHASED' predictor may be operated in two modes: either utilizing only 
the first phase and thereby producing two alternative dicer-cut location predictions, or utihzing 
both phases and thereby producing only one final dicer-cut location. 

In still another preferred embodiment of the present invention, the DICER-CUT 
LOCATION DETECTOR 116 preferably uses a 'FIRST-K' predictor, which utilizes the K- 
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nearest-neighbor algorifbm. The similarity metric between any two sequences is 1- E/L , where L 
is a parameter, preferably 8-10 and E is the edit distance between the two sequences, taking into 
account only the frst L nucleotides of each sequence. If the K-nearest-neighbor scores of two or 
more locations on the hairpin are not significantly different, these locations are further ranked by 
a Bayesian model, similar to the one described hereinabove. 

Scores of two or more of the abovementioned classifiers or predictors are 
integrated, yielding an integrated score for each 'potential VGAM'. As an example, Fig. 13C 
illustrates integration of scores from two classifiers, a 3' end recognition classifier and a 5' end 
recognition classifier, the scores of which are integrated to yield an integrated score. In a 
preferred embodiment of the present invention, INTEGRATED SCORE of 13C preferably 
implements a 'best-of-fareed' approach, accepting only 'potential VGAMs' that score highly on 
one of the above mentioned EDIT DISTANCE', or 'TWO-PHASED' predictors. Jn this context, 
'high scores' means scores which have been demonstrated to have low false positive value when 
scoring known miRNAs. 

The INTEGRATED SCORE is then evaluated as follows: (a) the 'potential 
VGAM' having the highest score is taken to be the most probable VGAM, and (b) if the 
integrated score of this 'potential VGAM' is higher than a pre-defined threshold, then the 
potential VGAM is accepted as the PREDICTED VGAM. 

Reference is now made to Fig. 7A which is a simplified block diagram of a 
preferred implementation of the TARGET-GENE BINDDING-SITE DETECTOR 118 described 
hereinabove with reference to Fig. 2. The goal of the TARGET-GENE BINDING-SITE 
DETECTOR 118 is to detect a BINDING SITE of Fig. 1, including binding sites located in 
untranslated regions of the RNA of a known gene, the nucleotide sequence of which BINDING 
SITE is a partial or accurate inversed reversed sequence to that of a VGAM RNA of Fig. 1, 
thereby determining that the above mentioned known gene is a host target gene of VGAM of Fig. 
1. 

The TARGET-GENE BINDING-SITE DETECTOR 118 of Fig. 2 therefore 
receives as input a plurality of DICER-CUT SEQUENCES FROM HAIRPINS 140 of Fig. 6A 
which were calculated by the previous step, and a plurality of POTENTIAL H OST TARGET 
GENE SEQUENCES 142 which derive from SEQUENCED DNA DATA 104 of Fig. 2, and after 
a phase of TARGET-GENE BINDING-SITE DETECTOR TRAINING & VALIDATION 128 is 
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operative to detect a plurality of POTENTIAL NOVEL TARGET-GENES HAVING BINDING 
SrrE/S 144 the nucleotide sequence of which is a partial or accurate inversed reversed sequence 
to that of each of the plurality of DICER-CUT SEQUENCES FROM HAIRPINS 140. Preferred 
operation of the TARGET-GENE BINDING-SITE DETECTOR 118 is further described 
hereinbelow with reference to Fig. 7B. 

Reference is now made to Fig, 7B which is a simplified flowchart illustrating a 
preferred operation of the target-gene binding-site detector 118 of Fig. 2. 

In a preferred embodiment of the present invention, the target-gene binding-site 
detector 118 first uses a sequence comparison algorithm such as BLAST in order to compare the 
nucleotide sequence of each of the plurality of dicer-cut sequences from hairpins 140, to the 
potential host target gene sequences 142, such a untranslated regions of known mRNAs, in order 
to find crude potential matches. Alternatively, the sequence comparison may preferably be 
performed using a sequence match search tool that is essentially a variant of the EDIT 
DISTANCE algorithm described hereinabove with reference to Fig. 6C, and the Nearest 
Neighbor algorithm (Mitchell, 1997). 

Results of tiie sequence comparison, perfomed by BLAST or other algorithms 
such as EDIT DISTANCE, are then filtered, preferably utilizing BLAST or EDIT DISTANCE 
score, to results which are similar to those of known binding sites (e.g. binding sites of miRNA 
genes Lui-4 and Let-7 to target genes Lin-14, Lin-41, Lin 28 etc.). Next the binding site is 
expanded, checking if nucleotide sequenced immediately adjacent to the binding site found by 
the sequence comparison algorithm (e,g, BLAST or EDIT DISTANCE), may improve the match. 
Suitable binding sites, then are confuted for free-energy and spatial structure. The results are 
analyzed, accepting only those binding sites, which have free-energy and spatial structure similar 
to that of known binding sites. Since known binding sites of known miRNA genes frequently 
. have multiple adjacent binding sites on the same target RNA, accordingly binding sites which are 
clustered are strongly preferred. Binding sites found in evolutionarily conserved sequences may 
preferably also be preferred. 

For each candidate binding site a score. Binding Site Prediction Accuracy, is 
calculated which estimates their similarity of its binding to that of known binding sites. This 
score is based on VGAM -binding site folding features including, but not limited to the free- 
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energy , the total number and distribution of base pairs, the total nmnb^ and distribution of 
unpaired nucleotides. 

la anotiier preferred embodiment of the present invention binding sites are 
searched by a reversed process: sequences of K(preferably 22)nucleotides of the untranslated 
regions of the target gene are assessed as potential binding sites. A sequence comparison 
algorithm, such as BLAST or EDIT DISTANCE, is then used to searcb. for partially or accurately 
complementary sequences elsewhere in the genome, which complementary sequences are found 
■ in known miRNA genes or computationally predicted VGAM genes. Only complementary 
sequences, the complementarity of which withstands the spatial structure and free energy analysis 
requirements described above are accepted. Clustered binding sites are strongly favored, as are 
potential binding sites and potential VGAM genes which occur in evolutionarily conserved 
genomic sequences. 

Host target binding sites, identified by the TARGET-GENE BINDING-SITE 
DETECTOR 118, are divided into 4 groiqjs: p) comprises binding sites that are exactly 
complementary to the predicted VGAM. a) b) and c) comprise binding sites that are not exactly 
complementary to the predicted VGAM: a) has binding sites with 0.9 < Binding Site Prediction 
Accuracy <= 1 ; b) has binding sites with 0.8 < Binding Site Prediction Accuracy <= 0.9 ; c) has 
binding sites witii 0.7 < Binding Site Prediction Accuracy <= 0.8 . The average number of 
mismatching nucleotides in the alignment of predicted VGAM and target binding site is smallest 
in category p) and largest in category c). 

In a preferred embodiment of the current invention a ranking of VGAM to host 
target gene binding is performed by calculating a score, Target Accuracy. This score is the 
dominant group identifier of all binding sites of a specific VGAM to a specific host target gene 
UTR, where 'a' dominates "b" and V dominates 'c'. 

In yet another preferred embodiment of the current invention a ranking of VGAM to host 
target gene binding is performed directly from the set of Binding Site Prediction Accuracies 
corresponding to all the binding sites of a specific VGAM to a specific host target gene UTR. 
This set of Accuracies is sorted in descending order. The final Target Accuracy is a sum of two 
terms: the first is a weighted sum of the sorted Accuracies where the weights are exponentially 
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decreasing as a function of the rank. The second team is a monotonously increasing function of 
the density of binding sites at the host target gene UTR. 

Host target binding genes, identified by the TARGET-GENE BINDING-STrB DETECTOR 118, 
are divided into 4 groups according to their target binding genes: A) 0.75 < Target Accuracy <= 
1 ; B) 0.65 < Target Accuracy <= 0.75 ; C) 0.5 < Target Accuracy <= 0.65 ; D) 0.3 < Target 

Accuracy <= 0.5 

Reference is now made to Fig. 8 which is a simplified flowchart illustrating a 
preferred operation of the function & utility analyzer 120 described hereinabove with reference to 
Fig. 2. The goal of the function & utility analyzer 120 is to detexmine if a potential host target 
gene is in fact a valid clinically useful target gene. Since a potential novel VGAM gene binding a 
binding site in the TJTR of a host target gene is understood to inhibit expression of that target 
gene, and if that host target gene is shown to have a valid clinical utility, then in such a case it 
follows that the potential novel viral gene itself also has a valid usefiil function - which is (he 
opposite of that of the host target gene. 

The function & utility analyzer 120 preferably receives as input a plurality of 
potential novel host target genes having binding-sit£^s 144, generated by the target-gene bmding- 
site detector 118, both of Fig. 7A. Each potential viral gene, is evaluated as follows: 

First, tiie system checks to see if the function of the potential host target gene is 
scientifically well established. Preferably, this can be achieved bioinformatically by searching 
various published data sources presentmg information on known function of proteins. Many such 
data sources exist and are published as is well known in the art. 

Next, for those host target genes the ftmction of which is scientifically known and 
is well documented, the system then checks if scientific research data exists which links them to 
known diseases. For example, a preferred embodiment of the present invention utiUzes the 
0M1M(TM) database published byNCBI, which summarizes research publications relating to 
genes which have been shown to be associated with diseases. 

Finally, the specific possible utility of the host target gene is evaluated. WhUe this 
process too may be facilitated by bioinformatic means, it mi^t require manual evaluation of 
published scientific research regarding the host target gene, in order to determine the utility of the 
host target gene to the diagnosis and or treatment of specific disease. Only potential novel viral 
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genes, the host target-genes of which have passed all three examinations, are accepted as novel 
viral genes. 



Reference is now made to Fig. 9, which is a simpUfied diagram describing each of 
a plurality of novel bioinformatically detected regulatory genes, referred to here as Viral 
Genomic Record (VGR) genes, which encodes an 'operon-Hke' cluster of novel micro RNA-like 
viral genes, each of which in turn modulates expression of at least one host target gene, the 
function and utility of which at least one host target gene is known in the art. 

VGR GE3!^ is a novel bioinformatically detected regulatory, non protein coding, 
RNA viral gene. The method by which VGR GENE was detected is described hereinabove with 
reference to Figs. 1-9. 

VGR GENE encodes VGR PRECURSOR RNA, an RNA molecule, typically 
several hundred nucleotides long. 

VGR PRECURSOR RNA folds spatially, forming VGR FOLDED PRECURSOR 
RNA. It is appreciated that VGR FOLDED PRECURSOR RNA comprises a plurality of what is 
known in the art as Tiaiipin' structures. These 'hairpin' structures are due to the fact that the 
nucleotide seque^ce of VGR PRECURSOR RNA comprises a plurality of segments, the first half 
of each such segment having a nucleotide sequence which is at least a partial or accurate inversed 
reversed sequence of the second half thereof, as is well known in the art. 

VGR FOLDED PRECURSOR RNA is naturally processed by ceUular enzymatic 
activity into separate VGAM precursor RNAs, herein schematically represented by VGAMl 
FOLDED PRECURSOR RNA through VGAM3 FOLDED PRECURSOR RNA, each of which 
VGAM precursor RNAs being a hairpin sh^ed RNA segment, corresponding to VGAM 
FOLDED PRECURSOR RNA of Fig.l . 

The above mentioned VGAM precursor RNAs are diced by DICER COMPLEX 
of Fig.l, yielding respective short RNA segments of about 22 nucleotides in length, 
schematically represraited by VGAMl RNA through VGAM3 RNA, each of which VGAM 
RNAs corresponding to VGAM RNA of Fig. 1 . 

VGAMl RNA, VGAM2 RNA and VGAM3 RNA, each bind complementarily to binding sites 
located in untranslated regions of respective target genes, designated VGAMl-TARGET RNA, 
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VGAM2-TARGET UNA and VGAM3-TARGET RNA, respectively, which target binding site 
conesponds to a target binding site such as BINDING SITE I, BINDING SITE n or BINDING 
SITE in of Fig. 1. This binding inhibits translation of the respective target proteins designated 
VGAMl-TARGET PROTEIN, VGAM2-TARGEr PROTEIN and VGAM3-TARGET 
PROTEIN respectively. 

It is appreciated that specific functions, and accordingly utiHties, of each VGR 
GENE of the present invention, correlates with, and may be deduced fiom, the identity of the 
target genes, which are inhibited by VGAM RNAs comprised in the 'operon-Hke' cluster of said 
VGR GENE, schematically represented by VGAMl TARGET PROTEIN through VGAM3 
TARGET PROTEIN. 

A listing of VGAM GENEs comprised in each of a plurality of VGR GENEs of 
Fig. 9 is provided in Table 7, hereby incorporated by refisience. Nucleotide sequ^ces of each 
said GAM GENEs and their respective genomic s ource and chromosomal 1 ocation are further 
described hereinbelow with reference to Table 1, hereby incorporated by reference. TARGET 
GENEs of each of said GAM GENEs are elaborated hereinbelow with reference to Table 4, 
hereby incorporated by reference. The functions of each of said TARGET GENEs and their 
association with various diseases, and accordingly the utilities of said each of GAM GENEs, and 
hence the functions and utiUties of each of said VGR GENEs of Fig. 9 is elaborated hereinbelow 
with reference to Table 5, hereby incorporated by reference. Studies establishing known 
functions of each of said TARGET GENEs, and correlation of each of said TARGET GENEs to 
known diseases are listed in Table 6, and are hereby incorporated by reference. 

The present invention discloses 205 novel genes of the VGR group of genes, 
which have been detected bioinformatically, as elaborated hereinbelow with reference to Table 7. 
Laboratory confirmation of 2 genes of the GR group of genes is described hereinbelow with 
reference to Figs. 14 through 17. 

In summary, the current invention discloses a very large number of novel VGR 
genes, each of which encodes a plurality of VGAM genes, whadi in turn may modulate 
expression of a plurality of host target proteins. 

Reference is now made to Fig. 10 which is a block diagram illustrating different 
utilities of genes of the novel group of genes of the present invention referred to here as VGAM 
genes and VGR genes. 
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The p resent i nvention discloses a first p lurality of n ovel v iral geaes r eferied t o 
here as VGAM geaes, and a second plurality of operon-like genes referred to here as VGR genes, 
eachoftheVGR genes encoding a plurality of VGAM genes. The present invention further 
discloses a very large number of known host target-genes, which are bound by, and the 
expression of which is modulated by each of the novel viial genes of the present invention. 
Published scientific data referenced by the present invention provides specific, substantial, and 
credible evidence that fhe above mentioned host target genes modulated by novel viral genes of 
the present invetition, are associated with various diseases. Specific novel genes of the present 
invention, host target genes thereof and diseases associated therewith, are described hereinbelow 
with reference to Tables 1 through 7. It is therefore appreciated that a function of VGAM genes 
and VGR genes of the present hivaation is modulation of expression of host target genes related 
to known diseases, and that therefore utilities of novel genes of the present invention include 
diagnosis and treatment of the above mentioned diseases. Fig. 10 describes various types of 
diagnostic and therapeutic utilities of novel genes of flie jwesent invention. 

A utility of novel genes of the present invention is detection of VGAM genes and 
of VGR genes. It is appreciated that since VGAM genes and VGR genes modulate expression of 
disease related host target g^aoes, that detection of expression of VGAM genes in clinical 
scenarios associated with said diseases is a specific, substantial and credible utility. Diagnosis of 
novel viral genes of the present invention may preferably be implemented by RNA expression 
detection tecbniques, including but not limited to biochips, as is well known in the art. Diagnosis 
of expression of genes of the present invention may be usefiil for research purposes, in order to 
fiirther understand the connection between the novel genes of the present invention and the above 
mentioned related diseases, for disease diagnosis and prevention purposes, and for monitoring 
disease progress. 

Another utility of novel viral genes of the present invention is and- VGAM gene 
therapy, a mode of therapy which allows up regulation of a disease related host target-gene of a 
novel VGAM gene of the present invention, by lowering levels of the novel VGAM gene which 
naturally inhibits expression of that host target gene. This mode of therapy is particularly useful 
with respect to target geaes which have been shown to be under-expressed in association with a 
specific d isease. A nti-VGAM gene t herapy i s further d iscussed h ereihbetow w ith r eference t o 
Figs. 1 lA and UB. 



31 



wo 2004/048511 PCT/IL2003/000998 

Reference is now made to Figs. IIA and IIB, simplified diagrams which when 
taken together illustrate anti-VGAM gene therapy mentioned hereinabove with reference to Fig. 
10. A utility of novel genes of the present invention is anti-VGAM gene ihecspy, a mode of 
ther^y which allows up regulation of a disease related host target-gene of a novel VGAM gene 
of the present invention, by lowering levels of the novel VGAM gene which naturally inhibits 
expression of that host target gene. Fig. 1 1 A shows a normal VGAM gene, inhibiting translation 
of a host target gene of VGAM gene, by binding to a BINDING SITE found in an untranslated 
region of VGAM TARGET RNA, as described hereinabove with reference to Fig. 1 . 

Fig. IIB shows an example of anti-VGAM gene therapy. ANTI-VGAM RNA is 
short artificial RNA molecule the sequence of which, is an anti-sense of V GAM RNA. Anti- 
VGAM treatment comprises transfecting diseased cells with ANTI-VGAM RNA, or with a DNA 
encoding thereof The ANTI-VGAM RNA binds the natural VGAM RNA, thereby preventing 
binding of natural VGAM RNA to its BINDING SITE. This prevents natural translation 
inhibition of VGAM TARGET RNA by VGAM RNA, thereby up regulating expression of 
VGAM TARGET PROTEIN. 

It is appreciated that anti-VGAM gene therapy is particularly usefiil with respect 
to host target genes which have been shown to be under-expressed in association with a specific 
disease. 

Furthermore, anti-VGAM therapy is particularly useful, since it may be used in 
situations in wMch technologies known in the art as RNAi and siRNA can not be utilized. As in 
known in the art, RNAi and siRNA are technologies which offer means for artificially inhibiting 
expression of a target protein, by artificially designed short RNA segments which bind 
complemientarily to mRNA of said target protein. However, RNAi and siRNA can not be used to 
directly upregulate translation of target proteins. 

Reference is now made to Fig. 12, which is a table summarizing laboratory 
validation results that validate efScacy of the bioinformatic gene detection engine 100 of Fig. 2. 
In order to assess efficacy of the bioinformatic gene detection engine 100, novel genes predicted 
thereby are preferably divided into 4 DETECTION ACCURACY GROUPS (first column), 
designated A through D, ranking VGAMS fi^om the most probable VGAMs to the least probable 
VGAMs, using the scores of HAIRPIN DETECTOR 114 and DICER-CUT LOCATION 
DETECTOR 116 as follows: 
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Groiq) A: The score of the HAIRPIN-DETECTOR is above 0.7, the overall score 
of the two-pliased predictor is above 0.55, and the score of the second phase of the two-phased 
predictor is above 0.75., or the score of the EDrT-DISTANCE predictor is equal or above 17. In 
this group, one miRNA is predicted for each haiipin. Group B: The score of the HAIRPIN- 
DETECTOR is above 0.5, the overall score of the two-phased predictor is above 0.55, and the 
hairpin is not in group A. Group C: The score of the HAIRPIN-DETECTOR is between 0.4 and 
0.5, and the overall score of the two-phased predictor is above 0.55. Group D: The score of the 
HAIRPIN-DETECTOR is between 0.3 and 0.4, and the overall score of the two-phased predictor 
is above 0.55. In groups B, C and D, if the score of the second phase of the two-phased predictor 
is above 0.75, one miRNA is predicted for each hairpin, otherwise both sides of the double 
stranded window are given as output, and are examined in the lab or used for binding site search. 
The groups are mutually exclusive, i.e. in groups A, C and D all hairpins score less than 17 in the 
EDir-DISTANCE predictor. 

It is appreciated fb.at the division into groups is not exhaustive: 410 of the 440 published 
hairpins (second column), and 896 of the 1560 novel VGAMs, belong to one of the groups. An 
indication o f t he r eal p erformance o f t he t wo p based p redictor i n t he p resence o f b ackground 
haiipii]^ is given by the column "precision on hairpin mixture' (third column). The precision on 
haiipin mixture is computed by mixing the published hairpins with background hairpins in a ratio 
of 1 :4 and taking as a woikmg assumption that they are hairpins not carrying a miRNA. This is a 
strict assun^tion, since some of these background hairpins may indeed contain iniRNAs, while 
in this column they are all coimted as failures 

Sample novel bioinformatically predicted human genes, of each of these groups 
are sent to the laboratory for validation (fourth column), and the number (fifth column) and 
percent (sixth column) of successful vaHdation of predicted human GAM genes is noted for each 
of the groups, as well as overall (bottom line). The number of novel VGAM genes exphcitly 
specified by present invention belonging to each of the four groups is noted (seventh column). 

It is ^predated that the present invention coni^rises 896 novel VGAM genes, 
which fall into one of these four detection accuracy groups, and that the bioinfoimatic gene 
detection engme 100 of Fig. 2 is substantiated by a group of 52 novel human GAM genes 
validated by laboratory means, out of 168 human GAM genes which were tested in the lab, 
resulting m validation of an overall 31% accuracy. The top group demonstrated 37% accuracy. 
Pictures of test-results of specific human genes in the abovementioned four groups, as well as the 
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mefliodology used for validating the expression of predicted genes is elaborated hereinbelow with 
refereacetoFig. 13. 

It is further appreciated that &ilure to detect a gene in the lab does not necessarily 
indicate a mistaken bioinfbnnatic prediction. Rather, it may be due to technical sensitivity 
limitation of the lab test, or because the gene is not expressed in fhe tissue examined, or at the 

development phase tested. 

It is still further appreciated that in general these findings are in agreement with 
the expected bioinformatic accuracy, as describe hereinabove with reference to Fig. 6B: assuming 
80% accuracy of fhe hairpin detector 114 and 80% accuracy of the dicer-cut location detector 1 16 
and 80% accuracy of the lab validation, this would result in 50% overall accuracy of the genes 
validated in the lab. 

Reference is now made to Fig. 13 which is a picture of labor atoiy results 
validating the eicpression of 37 novel human genes detected by the bioinformatic gene detection 
engine 100, in the four detection accuracy groups A through D described hereinabove with 
reference to Fig. 12. 

Each row in Fig. 13, designated A through D, correlates to a corresponding one of 
the four detection accuracy groups A-D, described hereinabove with reference to Fig. 12. Jn each 
row, pictures of several genes validated by hybridization of FCR-product southern-blots, are 
provided, each corresponding to a specific GAM gene, as elaborated hereinbelow. These PCR- 
product hybridization pictures are designated 1 through 22 in the A group, 1 throu^ 13 in the B 
group, 1 in the C group, and 1 in the D group. In each PCR hybridization picture, 2 lanes are 
seen: the test lane, designated and the control lane, designated For convenience of 
viewing the results, all PCR-product hybridization pictures of Fig. 13 have been shrunk x4 
vertically. It is ^preciated that for each of the tested genes, a clear hybridization band appears m 
the test ('+') lane, but not in the control ('-') lane. 

Specifically, Fig. 13 shows pictures of PCR-product hybridization validation by 
southem-blot, fhe methodology of which is described hereinbelow, to the following novel human 
GAM genes (RosettaGenomics Ltd. Gene Nomenclature): 

DETECTION ACCURACY GROUP A: (1) GAM8297.1; (2) GAM5346.1; 
(3)GAM281.1; (4) GAM8554.1; C5)GAM2071.1.; (6) GAM7553.1; (7) GAM5385.1; 
(8)GAM5227.1; (9) GAM7809.1; (10) GAM1032.1; (11) GAM3431.1; (12) QAM7933.1; 
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(13) GAM3298.1.;(14) GAM116.1; (15) GAM3418.1 (later published by other researchers as 
MIR23); (16) GAM3499.1; (17) GAM3027.1; (18) GAM7080.1; (19) GAM895.1; and (20) 
GAM2608.1, (21) GAM20, and (22) GAM21. 

DETECTION ACCURACY GROUP B: (1) GAM3770.1; (2) GAM1338.1; 
(3) GAM7957.1; (4) GAM391.1; (5) GAM 8678.1; (6) GAM2033.1; (7)GAM7776.1; (8) 
GAM8145.1; (9) GAM 633.1; (10) GAM19; (11) GAM8358.1; (12) GAM3229.1; and 
(13)GAM7052.1. 

DETECTION ACCURACY GROUP C: GAM25. 
DETECTION ACCURACY GROUP D:GAM7352.1. 
In addition to the PGR detection, the foUowiag GAMs were cloned and 
sequenced: GAM1338.1, GAM7809.1, GAM116.1, GAM3418.1 (later published by othea: 
researchers as MIR23), GAM3499.1, GAM3027.1, GAM7080.1, and GAM21. 

The PCR-pioduct hybridizatioii validation mefhodology used is briefly described 
as follows. In order to validate the expression of predicted novel GAMA^GAM genes, and 
assuming that these novel genes are probably expressed at low concentrations, a PGR product 
cloning approach was set up through the following strategy: two types of cDNA libraries 
designated "One tailed" and 'TLigation" were prepared fiom frozen HeLa SlOO extract (4c 
Biotech, Belgium) size fractionated RNA. Essentially, Total SlOO RNA was prepared through an 
SDS-Proteinase K incubation followed by an acid Phenol-Chloroform purification and 
Isopropanol precipitation. Alternatively, total HeLa RNA was also used as starting material for 
these libraries. 

Fractionation was done by loading up to 500ng per YMIOO Amicon Microcon 
column (MiUipore) followed by a 500g centrifiagatioii for 40 minutes at 4"C. Flowthrough 
"YMIOO" RNA consisting of about % of the total RNA was used for library preparation or 
fractionated further by loading onto a YM30 Amicon Microcon column (Milhpore) followed by 
a 13,500g centrifugation for 25 minutes at 4''C, Flowthrough "YM30" was used for library 
preparation as is and consists of less than 0.5% of total RNA. For the both the "ligation" and the 
"One-tailed" libraries RNA was dephosphorilated and ligated to an RNA (lowercase)-DNA 
(UPPERCASE) hybrid 5'-phosphoiilated, 3' idT blocked 3 '-adapter (5'-P- 
uuuAACCGCATTCTC-idT-3' Dhannacon # P-002045-01-05) (as elaborated in Elbashir et al 
2001) resulting in ligation only of RNase m type cleavage products. 3'-Ligated RNA was 
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excised and purified firom a half 6%, half 13% polyacrylamide gel to remove excess adapter with 
a Nanosep 0.2|iM centrifiigal device (Pall) according to instructions, and precipitated with 
glycogen and 3 volumes of Ethanol. Pellet was resuspended in a minimal volume of water. 

For the "ligation" library a DNA (UPPERCASE)-RNA (lowercase) hybrid 5'-adapter (5'- 
TACTAATACGACTCACTaaa-3' Dhannacon # P-002046-01-05) was ligated to the 3'-adapted 
FNA, reverse transcribed with "EcoRI-RT" (5'- 

GACTAGCTGGAATTCAAGGATGCGGTTAAA-3'), PGR amplified with two external 
primers essentially as in Elbashir et al 2001 except that primers were "EcoRI-RT" and 'Tsti 
Fwd" (5'-CAGCCAACGCTGCAGATACGACTCACTAAA-3'). This PGR product was used as 
a template for a second round of PGR with one hemispecific and one external primer or with two 
hemispeoific primers. 

For the "One tailed" library the 3 '-Adapted RNA was annealed to 20pmol primer 
"EcoRI RT' by heating to IGTC and cooling 0.1°C/sec to 30^ and then reverse transcribed with 
Superscript n RT (According to instructions, Invitrogen) in a 20jil volume for 10 alternating 5 
minute cycles of 3TC and 45''C. Subsequently, KNA was digested. with lul 2M NaOH, 2mM 
EDTA at eS'C for 10 minutes. cDNA was loaded on a polyacrylamide gel, excised and gel- 
purified firom excess primer as above (invisible, judged by primer run alongside) and 
resuspended in 13(il of water. Purified cDNA was then oligo-dC tailed with 400U of recombinant 
terminal transferase (Roche molecular biochemicals), Ifil lOO^M dCTP, l]xl 15mM C0CI2, and 
4nl reaction buffer, to a final volume of 20nl for 15 minutes at SVC. Reaction was stopped with 
2jil 0.2M EDTA and 15|il 3M NaOAc pH 52. Volume was adjusted to 150jal with water, Phenol 
: Bromochloropropane 10:1 extracted and subsequently precipitated with glycogen and 3 
volumes of Ethanol. C-tailed cDNA was used as a template for PGR with the external primers 
'T3-PsfBsg(G/I)i8" 

(5'-AATTAACCCTCACTAAAGGCTGCAGGTGCAGGIGGGnGGGnGGGnGN-3' where I 
stands for Ihosine and N for any of the 4 possible deoxynucleotides), and with "EcoRI Nested" 
(5'-GGAATTCAAGGATGCGGTTA-3'). This PGR product was used as a template for a second 
round of PGR with one hemispecific and one external primer or with two hemispecific primers. 

Hemispecific primers were constructed for each predicted GAM/VGAM by an in- 
house program designed to choose about half of the 5' or 3' sequence of the GAftWGAM 
corresponding to a TM" of about 30''-34°C constrained by an optimized 3' clamp, impended to the 
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cloning adaptec s equence (for " One-tailed" libraries 5 '-GGMNGGGNNG on the 5 ' end o f the 
GAMA^GAM, or TTTAAGCGCATC-3' on the 3' end of the GAMATGAM. For "Ligation" 
libraries the same 3' adapter and 5'-CGACTCACTAAA on the 5' end). Consequently, a fully 
complementary primer of a TM" higher than 60°C was created covering only one half of the 
GAMA^GAM sequence permitting the unbiased elucidation by sequencing of the other half. 

CONFIRMATION OF GAMA^GAM SEQUENCE AUTHENnCITY OF PCRPRODUCTS: 

SOUTHERN BLOT: PCR-prcduct sequences were confirmed by southern blot 
(Soutiiem EM. Biotechnology.l992;24: 122-39.(1975)) and hybridization with DNA 
oligonucleotide probes synthesized against predicted GAMsA^GAMs. Gels were transierred onto 
a Biodyne PLUS 0.45 [im, (PaH) positively charged nylon membrane and UV cross-linked. 
Hjiiridization was performed ovemigjit with DIG-labeled probes at 42°C in DIG Easy-Hyb 
buffer (Roche). Membranes were washed twice with 2xSSC and 0.1% SDS for 10 min. at 42''C 
and then washed twice with O.SxSSC and 0.1% SDS for 5 min at 42°C. The membrane was then 
developed by using a DIG lumiaescent detection kit (Roche) using anti-DIG and CSPD reacdon, 
according to the manufecturer's protocol. All probes were prepared according to the 
manufacturers (Roche Molecular Biochemicals) protocols: Digoxigenin (DIG) labeled antisense 
transcripts was prepared from purified PGR products using a DIG RNA labeling kit with T3 
RNA polymerase. DIG labeled PGR was prepared by using a DIG PGR labeling kit. 3'-DIG- 
tailed oligo ssDNA antisense probes, containing DIG-dUTP and dATP at an average tail length 
of 50 nucleotides were prepared &om lOOpmole oligonucleotides with the DIG Oligonucleotide 
Labeling Kit. 

CLONING: PGR products were inserted into pGEM-T (Promega) or pTZ57 
(MBI Fermentas), transformed into competent JM109 £. coli (Promega) and sown on LB-An^ 
plates with IPTG/Xgal. White and light-blue colonies were transferred to duplicate gridded 
plates, one of which was blotted, onto a membrane (Biodyne Plus, Pall) for hybridization with 
DIG tailed oligo probes (according to instructions, Roche) corresponding to the expected GAM. 
Plasmid DNA from positive colonies was sequenced. 

Reference is now made to Fig. 14A, which is a schematic representation of a novel 
human GR gene, herein designated GR12731 (RosettaGenomics Ltd. Gene Nomenclature), 
located on chromosome 9, comprising 2 known MIR genes - MIR23 MIR24, and 2 novel GAM 
genes, herein designated GAM22 and GAMl 16, all marked by solid black boxes. Fig. 14A also 
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schematically illustrates 6 non-GAM hairpin sequences, and one non-haitpin sequence, all 
marked by white boxes, and serving as negative controls. By 'non-GAM hairpin sequences* is 
meant sequences of a similar length to known MIR PRECURSOR sequences, which form hairpin 
secondary folding pattern similar to MIR PRECURSOR hairpms, and yet which are assessed by 
the biomfonnatic gene detection engine 100 not to be valid GAM PRECURSOR haiipms. It is 
appreciated that Fig. 14A is a sunpli£ed schematic representation, reflecting only the order in 
which the segments of interest appear relative to one another, and not a proportional distance 
between the segments. 

Reference is now made to Fig. 14B, which is a schematic representation of 
secondary folding of each of the MIRs and GAMs of GR GR12731 - MIR24, MIR23, GAM22 
and GAM116, and of the negative control non-GAM hairpins, herein designated N2, N3, N116, 
N4, N6 and N7. NO is a non-hairpin control, of a sunilar length to that of known MIR 
PRECURSOR hairpins. It is appreciated that the negative controls are situated adjacent to and in 
between real MIR and GAM genes, and demonstrates sunilar secondary folding pattcms to that 
of known MIRs and GAMs. 

Reference is now made to Fig. 14C, which is a picture of laboratory results of a 
PGR test upon a YMIOO "ligation"-librafly, utilizing specdjSc primer sets directly inside the 
boundaries of the hairpins. Due to the nature of the library the only PGR amphfiable products 
can result fiom RNaselEI type enzyme cleaved RNA as expected for legitimate hairpin 
precursors presumed to be produced by DROSHA (Lee et al, Nature 425 415-419, 2003). Fig 
14C demonstrates expression of hairpin precursors of known MIR genes - MIR23 and MIR24, 
and of novel bioinfonnatically detected GAM22 and GAM116 genes predicted bioinformatically 
by a system constructed and operative in accordance with a preferred embodiment of the present 
invention. Fig. 14C also shows that none of the 7 controls (6 hairpins designated N2, N3, N23, 
N4, N6 and N7 and 1 non-hairpin sequence designated NO) were expressed. Nl 16 is a negative 
control sequence partially overlapping GAM116. 

In the picture, test lanes including template are designated and the control lane is designated 
It is appreciated that for each of the tested hairpins, a clear PCR band appears in the test ('+') 
lane, but not in the control ('-') lane. 

Figs. 14A through 14C, when taken together validate the efificacy of the 
bioinformatic gene detection engine in: (a) detecting known MIR genes; (b) detectmg novel 
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GAM genes which are found adjacent to these MIR genes, and which despite exhaustive prior 
biological efforts andbioinfonnatic detection effi)rts, went undetected; (c) disceining between 
GAM (or MIR) PRECURSOR hairpins, and non-GAM hairpins. 

It is jrppreciated that the ability to discern GAM-haiipins from non-GAM-hairpins 
is very significant in detecting GAM genes, since hairpins in general are highly abundant in ■flie 
genome. Other MIR prediction programs have not been able to address this challenge 
successfiilly. 

Reference is now made to Fig. 15A which is an annotated sequence of human EST 
comprising a novel gene detected by the gene detection system of the present invention. Fig. 
15A shows the nucleotide sequence of a known human non-protein coding EST (Expressed 
Sequence Tag), identified as EST72223. The EST72223 clone obtained fi:Dm TIGR database 
(Kirkness and Kerlavage, 1997) was sequenced to yield the above 705bp transcript with a 
polyadenyl taU. It is ^predated that the sequence of this EST comprises sequences of one known 
nuRNA gena, identified as MIR98, and of one novel human GAM gene, refetred to here as 
GAM25, detected by the bioinfoimatlc gene detection system of the present invention and 
described hereinabove with reference to Fig. 2. 

The sequences of the precursors of the known MIR98 and of the predicted 
GAM25 are in bold, the sequences of the established miRNfA 98 and of the predicted miRNA 
GAM25 are underlined. 

Reference is now made to Figs. ISB, 15C and 15D that are pictures of laboratory 
results, which when taken together demonstrate laboratory confirmation of expression of the 
bioinformatically detected novel gene of Fig. 15A. 

In two parallel experiments, an enzymatically synthesized capped, EST72223 
RNA transcript, was incubated with Hela 8100 lysate for 0 minutes, 4 hours and 24 hours. RNA 
was subsequently harvested, run on a denaturing polyacrylamide gel, and reacted with a 102nt 
and a I45nt antisense MIR98 and GAM25 precursor transcript probes respectively. The Northem 
blot results of these experiments demonstrated processing of EST72223 RNA by Hela lysate 
(lanes 2-4, in 15B and 15C), into ~80bp and ~22bp segments, which reacted with the MIR98 
precursor probe (15B), and into -lOObp and ~24bp segments, which reacted with the GAM25 
precursor probe (15C). These results demonstrate the processing of EST72223 by Hela lysate into 
MIR98 precursor and GAM25 precursor. It is also appreciated firom Fig. 15C (lane 1) that Hela 
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lysate itself reacted with the GAM25 precursor probe, in a number of bands, including a ~100bp 
band, indicating that GAM25-precursor is endogenously expressed in Hela cells. The presence of 
additional bands, higher than lOObp in lanes 5-9 probably corresponds to the presence of 
nucleotide sequences in Hela lysate, which contain the GAM25 sequence. 

In addition, in order to demonstrate the kinetics and specificity of the processing 
of MIR98 and GAM25 nuRNA precursors into their respective miRNA's, transcripts of MIR98 
and of the bioinformatically predicted GAM25, were sknilarly incubated with Hela SI 00 lysate, 
for 0 mitiutes, 30 minutes, 1 hour and 24 hours, and for 24 hours with the addition of EDTA, 
added to inhibit Dicer activity, following which RNA was harvested, run on a polyacrylamide gel 
and reacted with MIR98 and GAM25 precursor probes. Capped transcripts were prepared for in- 
vitro RNA cleavage assays with T7 RNA polymerase including a m'G(5')ppp(5')G-cappiQg 
reaction using the mMessage luMachine kit (Ambion). Purified PGR products were used as 
template for the reaction. These were amplified for each assay with specific primers containing a 
T7 promoter at the 5' end and a T3 RNA polymerase promoter at the 3'end. Capped RNA 
transcripts were incubated at 30°C in supplemented, dialysis concentrated, Hela SlOO 
cytoplasmic extract (4C Biotech, Seneffe, Belgium). The Hela SlOO was supplemented by 
dialysis to a final coucentratioD of 20mM Hepes, lOOmM KCl, 2.SmM MgCk , O.SmM DTT, 
20% glycerol and protease inhibitor cocktail tablets (Complete mini Roche Molecular 
Biochemicals). After addition of all components, final concentrations were lOOmM capped 
target RNA, 2mM ATP, 0.2mM GTP, 500U/ml RNasin, 25p.g/ral creatine kinase, 25mM creatine 
phosphate, 2.5mM DTT and 50% SlOO extract. Proteinase K, used to enhance Dicer activity 
(Zhang H, Kolb FA, Brondani V, Billy E, Filipowicz W. Human Dicer preferentially cleaves 
dsRNAs at their tennini without a requirement for ATP.EMBO J. 2002 Nov l;21C21):5875-85) 
was dissolved in SQmM Tris-HCl pH 8, SniM CaCl2, and 50% glycerol, was added to a fibial 
concentration of 0.6 mg/ml. Cleavage reactions were stopped by the addition of 8 volumes of 
proteinase K buffer (200Mm Tris-Hcl, pH 7.5, 25m M EDTA, 300mM NaCI, and 2% SDS) and 
incubated at 65°C for 15niin at different time points (0, 0.5, 1, 4, 24h) and subjected to 
phenol/chloroform extraction. Pellets were dissolved in water and kept frozen. Samples were 
analyzed on a segnented half 6%, half 13% polyacrylamide 1XTBE-7M Urea gel. 

The Northern blot results of these experiments demonstrated an accumulation of a 
~22bp segment which reacted with the MIR98 precursor probe, and of a ~24bp segment which 
reacted with the GAM25 precursor probe, over time (lanes 5-8). Absence of these segments when 
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incubated with EDTA (lane 9), which is known to inhibit Dicer eaizyme (Zhang et al., 2002), 
supports the notion that the processing of MERdS and GAM25 miRNA's firom their precursors is 
mediated by Dicer enzyme, found in Hela lysate. The molecular sizes of BST72223, MIR-98 and 
GAM25 and their corresponding precursors are indicated by arrows. 

Fig. 15D present Northern blot results of same above experiments with GAM25 probe 
(24nt). The results clearly demonstrated the accumulation of mature GAM25 gene after 24 h. 

To validate the identity of the band shown by the lower arrow in Figs. 15C and 15D, a 
SNA band parallel to a marker of 24 base was excised firom the gel and cloned as in Elbasfair et al 
(2001) and sequenced. 90 clones corresponded to the sequence of mature GAJM25 gene, three 
corresponded to GAM25* (the opposite arm of the haiipin with a 1-3 nucleotide 3' overhang) 
and two to the hairpin-loop. 

GAM25 was also validated endogenously by sequencing from both sides from HeLa 
YMIOO total-RNA "ligation" hbraries, utilizing hemispecific primers as detailed in Figl3. 

Taken together, these results validate the presence and processing of a novel MIR 
gene product, GAM25, which was predicted bioinfonnatically. The processing of this novel gene 
product, by Hela lysate from EST72223, through its precursor, to its jSnal form was similar to 
that observed for known gene, MIR98. 

Transcript products were 705nt (EST72223), I02nt (MIR98 precursor), 125nt (GAM25 
precursor) long. EST72223 was PGR amplified with T7-EST 72223 forward primer: 

5*-TAATACGACTCACTATAGGCCCTTATTAGAGGATTCTGCT-3' 

and T3-EST72223 reverse primer: 

5'-AATTAACCCTCACTAAAGGTTTTTTTTTCCTGAGACAGAGT-3'. 
MIR98 was PGR amplified using EST72223 as a template with T7MIR98 forward primer: 
5-'TAATACGACTCACTATAGGGTGAGGTAGTAAGTTGTATTGTT-3' 
and T3MIR98 reverse primer: 



5'-AATTAACCCTCACTAAAGGGAAAGTAGTAAGTTGTATAGTT-3'. 
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GAM25 was PGR amplified using EST72223 as a template with GAM25 forward primer: 5'- 
GAGGCAGGAGAATTGCTTGA- 3' and T3-EST72223 reverse primer: 



5'-AATTAACCCTCACTAAAGGCCTGAGACAGAGTCTTGCTC-3'. 

It is appreciated that the data presented in Figs, 15A, 15B, 15C and 15D when 
taken together validate the fimction of the bioinformatic gene detection engine 100 of Fig. 2. Fig. 
15A shows a novel GAM gene bioinformatically detected by the bioinformatic gene detection 
engine 100, and Figs. 15C and 15D show laboratory confirmation of the expression of this novel 
gene. This is in accord with the engine training and validation methodology described 
hereinabove with reference to Fig. 3, 

It is appreciated by persons skilled in the art that the present invention is not 
limited by what has been particularly shown and described hereinabove. Rather the scope of the 
present invention includes both combinations and subcombinations of the various features 
described hereinabove as well as variations and modifications which would occur to persons ' 
skilled in the art upon reading the specifications and which are not in the prior art. 

Reference is now made to Fig. 16, which, presents pictures of laboratory results, 
which demonstrate laboratory confirmation of excision ("dicing") of the bioinformatically 
detected novel VGAM HTVl genes, herein designated VGAM2032.2, VGAM3249.1, VGAM 
507.2 and VGAM1016.2 firom their predicted precursors by their incubation in HeLa S-100 
lysate as described in Fig. 15. 

Fig.l6A presents the entire 5'-UTR of HTVl (USR) containing two predicted VGAM 
precursor genes, in bold; VGAM 2032 and VGAM3249. The bioinformatically predicted mature 
VGAMs are depicted with underscore, the 5 '-most is VGAM 2032.2 and the second is VGAM 
3249. 1 . VGAM 2032.2 matches the known HIVl RNA structure named TAR to which the TAT 
protein binds OTature 1987. 330:489-93). 

Figl6B and Fig 16C depict Northern blot analysis of VGAMs in U5R, hybridized with 
predicted mature VGAM oligonucleotide probes VGAM 2032.2, and VGAM 3249.1, 
respectivley. The molecular size of the entire USR transcript, 355nt, is indicated by arrow. The 
predicted molecular sizes of VGAM 2032.2, and VGAM 3249.1 axe 22nt and 17ht respectively. 
The 22nt molecular madter is indicated by arrow. Lanes: 1-Hela lysate; 2-U5R transcript in HeLa 
Lysate without incubation. 3- U5R transcript incubated overnight with Hela lysate. 
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Figs.l6D and 16E present partial transcripts of HIVl RNA reacted with predicted mature 
HIVl-VGAM oligoniicleotide probes. lii each figure, the experimentait transcript sequence is 
shown, annotated in bold is the predicted VGAM precursor, and in underscore the predicted 
mature VGAM. Norflieni blot analysis of VGAM precursors for VGAM507.2 (Figl6D), and 
VGAM1016.2 (Fig.l6E). The transcaipt sizes are 163nt for VGAM507.2 transcript and 200nt for 
VGAM1016.2 transcript. The predicted molecular sizes of VGAM507.2 and VGAM1016.2 are 
both 24nt The 22nt molecular marker is indicated by arrow. Lanes: 1- Transcript in HeLa Lysate 
without incubation. 2- Transcript incubated overnight with HeLa lysate. 

It is £5)preciated that the sequence of the reacting bands in the foreseen sizes comprise 
sequences of novel VGAM genes, referred to here as VGAM 2032.2, VGAM 3249.1, VGAM 
507.2 , VGAM 1016.2 , detected by the bioinformatic gene detection engine ICQ of the present 
invention, described herekiahove with reference to Fig. 2. 

Reference is now made to Fig. 17 which presents pictures of laboratory results, which 
demonstrate laboratory confiimation of expression of the bioinfonnatically detected novel 
Vaccinia VGAM genes VGAM224 (Figs. 17A and 17C) and VGAM3184 (Fig.l7B). HeLa cells 
were infected with 50 PFU Vaccinia Vmis and total RNA was harvested after 3 days. Northern 
blot analysis of VGAM precursors in total RNA extracted from HeLa cells infected with 
Vaccinia Virus, lane 1, or HeLa uninfected cells, lane 2, and hybridized with predicted precursor 
DIG-labled RNA probes of 53nt for VGAM224 (Fig. 17A), or of 73nt for VGAM3184 (Fig. 
17B) or with a 22nt y'^P-ATP-Iabled DNA oligo probe for predicted mature VGAM224.2 (Fig. 
17C). A transcript of predicted sequence and size was run alongside as a size marker and as a 
hybridization control, lane 3 (except Fig.l7C). Arrow in Fig. 17C marks band of expected 
precursor size, 53nt, reactmg with mature 22nt VGAM224.2 probe. 

It is appreciated that the sequence of the reacting bands appear only in infected cells in- 
vivo and comprise sequences of novel VGAM gene precursors, refened to here as VGAM224 
and VGAM3184, detected by tiie bioinformatic gene detection engine 100 of the present 
invention, as described hereinabove with reference to Fig. 2. 

DETAILED DESCRIPTION OF LARGE TABLES 

Table 1 comprises data relating to the source and location of novel VGAM genes 
of the present invention, and contains the following fields: 

GENE NAME Rosetta Genomics Ltd. gene nomenclature (see below) 
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VGAM SEQ-ID VGAM Seq-ID, as in the Sequence listing 

PRECUR SEQ-ID VGAM precursor Seq-ID, as in the Sequence Listing 

ORGANISM Virus name 

GENOME TYPE Genome type of the virus; dsRNA, dsDNA, ssRNA negative- 
stiaud, ssRNA positive-strand, Deltavirus or Retroid, as taken 
firom ORGANISM definition by GenBank, NCBI. 

GENOME STRUCTURE genome organization: circular or linear. 

SOURCE_REF-ID Accession number of vims source sequence 

SOURCE_OFFSET Offset of VGAM precuisor sequence on source sequence 

STRAND (+) positive strand, (-) negative strand 

SRC Source-type of VGAM precursor sequence (see below) 

VGAM ACC VGAM Prediction Accuracy Group (see below); 

Table 2 comprises data relating to VGAM precursors of novel VGAM genes of 
the present invention, and contains the following fields: 

GENE NAME Rosetta Genomics Ltd. gene nomenclature (see below) 

PRECUR SEQ-ID VGAM precursor Seq-ID, as in the Sequence Listing 

PRECURSOR SEQUENCE VGAM precursor nucleotide sequence (5' to 3') 

FOLDED-PRECURSOR Schematic representation of the VGAM folded precursor, 

beginning 5' end (beginning of upper row) to 3' end 
(beginning of lower row), where the hairpin loop is 
positioned at the right part of the draw. 

SRC Source-type of VGAM precursor sequence (see below) 

VGAM ACC VGAM Prediction Accuracy Group (see below); 

Table 3 comprises data relating to VGAM genes of the present invention, and 
contains the followmg fields: 

GENE NAME Rosetta Genomics Ltd. gene nomenclature (see below) 

VGAM SEQ-ID; VGAM Seq-ID, as in tiie Sequence Listing 
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GENE_SEQUENCE Sequence (5' to 3') of the mature, *diced' VGAM gene 

PRBCUR SEQ-ID VGAM precursor Seq-DD, as in the Sequence Listing 

SOURCE_KEF-n) Accession number of the source sequence 

SRC Source-type of VGAM precursor sequence (see below) 

VGAM ACC VGAM Prediction Accuracy Group (see below); 



Table 4 comprises data relating to host target-genes and binding sites of VGAM 
genes of the present invention, and contains the following fields: 

GENE NAME Rosetta Genomics Ltd. gene nomenclature (see below) 

VGAM SEQ-ID; VGAM Seq-ID, as in the Sequence Listing 

TARGET VGAM target protem name 

#BS Number of unique binding sites of VGAM onto Target 

TARGET SEQ-ID Target binding site Seq-ID, as in the Sequence Listing 

TARGET REF-ID Target accession number (GenBauk) 

JJTR Untranslated region of binding site/s (3 ' or 5 ') 

UTR OFFSET Offset of VGAM binding site relative to UTR 

TAR-BS-SEQ Nucleotide sequence (5' to 3') of the host target binding site 

BINDING-SITE-DRAW Schematic representation of the binding site, upper row 

present 5 ^ to 3 ' sequence of the VGAM, lower row present 3 ' 
to 5" sequence of the target. 

SRC Source-type of VGAM precursor sequence (see below) 

VGAM ACC VGAM Prediction Accuracy Group (see below); 

BS ACC Binding-Site Accuracy Group (see below) 

TAR ACC Target Accuracy Group (see below); 

Table 5 comprises data relating to functions and utilities of novel VGAM genes of 
the present invention, and contains the following fields: 
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GENE NAME Rosetta Genomics Ltd. gene nomenclature (see below) 

TARGET VGAM target protein name 

GENE_SEQUENCE Sequence (5' to 3') of the mature, 'diced' VGAM gene 

GENE-FtfNCnON Description of the VGAM functions and utilities 

SRC SouFce-type of VGAM precursor sequence (see below) 

VGAM ACC VGAM Prediction Accuracy Group (see below) 

TAR ACC Target Accuracy Group (see below) 

TAR DIS Target Disease Relation Group (see below) 

Table 6 comprises a bibliography of references supporting the fimctions and 
utilities of novel VGAM genes of the present invention, and contains the following fields: 
GENE NAME Rosetta Gnomics Ltd. gene nomenclature (see below) 

TARGET VGAM target protein name 

REFERENCES list ofreferences relating to the host target gene, 

SRC Source-type of VGAM precursor sequence (see below) 

VGAM ACC VGAM Prediction Accuracy Group (see below) 

TAR ACC Target Accuracy Group (see below); and 

Table 7 comprises data relating to novel VGR genes of the present invention, and 
contains the following fields: 

GENE NAME ' Rosetta Genomics Ltd. VGR gene nomenclature 

SOURCE START OFFSET Start-oflfeet of VGR gene relative to source sequence 
SOURCE END OFFSET End-offset of VGR gene relative to source sequence 
SOURCE_REF-ID Accession number of the source sequence 

STRAND (+) positive strand, (-) negative strand 

VGAMS_ID'S_IN_VGR List of the VGAM genes in the VGR. cluster 
SRC Source-type of VGAM precursor sequence (see below) 

VGR ACC VGR Prediction Accuracy Group (see below). 
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The following conventions and abbreviations are used in the tables: 

GENE NAME is a RosettaGeaoimcs Ltd. gene nomenclature. All VGAMs are 
designated by VGAMx. 1 or VGAMx.2 where x is the unique SEQ-ID. If the VGAM precursor 
has a single prediction for VGAM, it is designated by VGAMx. 1. Otherwise, the higher accuracy 
VGAM prediction is designated by VGAMx. I and the second is designated by VGAMx.2. 

SRC is a field indicating the type of source in which novel genes were detected, as 
one of the following options: (100) DNA sequence, (101) RNA sequence. Sequences are based 
on NCBI Build33 of the viral genome annotation. 

VGAM ACC (VGAM Prediction Accuracy Group) of gene prediction system: A- 
very high accuracy, B- high accuracy, C- moderate accuracy, D-low accuracy, as described 
hereinbelow with reference to Flg.l2. 

BS ACC (Binding-Site Accuracy Gxovp) indicates accuracy of target bmding site 
prediction. A- very high accuracy, B- high accuracy, C- moderate accuracy, as described 
hereiribelow with reference to Fig.l4B. 

TAR ACC (Target Accuracy Group) indicates accuracy of total GAM-target 
binding prediction, considering the number of bindiog sites a GAM has on the target's UTR; A- 
very high accuracy, B- hi^ accuracy, C- moderate accuracy, as described hereinbelow with 
reference to Fig.l4B. 

TAR DIS (Target Disease Relation Group) 'A' indicates if the target gene is 
known to have a specific causative relation to a specific known disease, based on the OMIM 
database. It is appreciated that this is a partial classification emphasizing genes which are 
associated with 'smgle gene' diseases etc. All genes of the present invention ARE associated 
with various diseases, although not all are in 'A' status. 

VGR ACC (GR Prediction Accuracy Group) indicates the maximum gene 
prediction accuracy among VGAM genes of the cluster. A- very high accuracy, B- high accuracy, 
C- moderate accuracy, as described herdhbelow with reference to Fig.l4B. 
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CLAIMS 

1 . A bioinfoimatically detectable novel viral gene encoding substantially pure 
nucleic acid wherein: 

RNA encoded by said bioinformatically detectable novel viral gene is about 18 to about 
24 nucleotides in length, and originates from an RNA precursor, which RNA precursor is about 
50 to about 120 nucleotides in length; 

a nucleotide sequence of a first half of said RNA precursor is a partial inversed-reversed 
sequence of a nucleotide sequence of a second half thereof; 

a nucleotide sequence of said RNA encoded by said novel viral gene is a partial inversed- 
reversed sequence of a nucleotide sequence of a binding site associated with at least one host 
target gene; and 

a function of said novel viral gene is bioinformatically deducible. 

2. A bioinformatically detectable novel viral gene encoding substantially pure DNA 
wherein: 

RNA encoded by said bioninformatically detectable novel viral gene comprises a plurality 
of RNA sections, each of said RNA sections being about 50 to about 120 nucleotides in length, 
and comprising an RNA segment, which RNA segment is about 18 to about 24 nucleotides in 
length; 

a nucleotide sequence of a first half of each of said RNA sections encoded by said novel 
gene is a partial inversed-reversed sequence of nucleotide sequence of a second half thereof; 

a nucleotide sequence of each of said RNA segments encoded by said novel gene is a 
partial inversed-reversed sequence of the nucleotide sequence of a binding site associated with at 
least one target gene; and 
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a function of said novel gene is bioinfonaatically deducible JEbom the following data 
elements: said nucleotide sequence of said RNA encoded by said novel gene, a nucleotide 
sequence of said at least one target gene, and fiinction of said at least one target gene. 

3 . A bioinfonnatically detectable novel viral gene encoding substantially pure DNA 
wherein: 

RNA encoded by said bioin&imatically detectable novel gene is about 18 to about 24 
nucleotides in length, and originates from an RNA precursor, which RNA precursor is about 50 
to about 120 nucleotides in length; 

a nucleotide sequence of a first half of said RNA precursor is a partial inversed-reversed 
sequence of a nucleotide sequence of a second half thereof; 

a nucleotide sequence of said RNA encoded by said novel gene is a partial inversed- 
reversed sequence of a nucleotide sequence of a binding site associated with at least one target 
gene; 

a function of said novel gene is modulatioa of expression of said at least one taiget gene; 

and 

said at least one target gene does not encode a protein. 

4. A bioinfoimatically detectable novel gene according to claim 1 and wherein said 
function of said novel gene is bicinfonnatically deducible from the following data elements: 

said nucleotide sequence of said RNA encoded by said bioinformatically detectable novel 

gene, 

a nucleotide sequence of said at least one target gene; and 

a fimction of said at least one target gene. 

5. A bioinformatically detectable novel gene according to claim 1 and wherein said 
RNA encoded by said novel gene complementarily binds said binding site associated with said at 
least one target gene, thereby modulating expression of said at least one target gene. 
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6. A bioinformaticaUy detectable novel gene according to claim 1 and wherein: 

said binding site associated with at least one target gene is located in an untranslated 
region of RNA encoded by said at least one taarget gene. 

7. A bioinftumatically detectable novel gene according to claim 5 and wherein: 

said fimction of said novel gene is selective inhibition of translation of said at least one 
target gene, which selective iohibition comprises complementary hybridization of said RNA 
encoded by said novel gene to said binding site. 

8. A vector comprising the DNA of claim 1 . 

9. A method of selectively inhibiting translation of at least one gene, comprising 
introducing the vector of claim 8 into a cell. 

10. A method according to claim 9 and wherein said introducing comprises utilizing 
RNAi patiiway. 

11. A gene expression inhibition system comprising: 
the vector of claim 8; and 

a vector inserter, functional to insert said vector of claim 10 into a cell, thereby selectively 
inhibiting translation of at least one gene. 



12. A probe comprising the DNA of claim 1 . 

13. A method of selectively detecting expression of at least one gene, comprising 
using the probe of claim 12(14). 

14. A gene expression detection system comprising: 
the probe of claun 12; and 

a geae expression detector functional to selectively detect expression of at least one gene. 

15. An anti-viral substance capable of neutralizing said RNA of claim 1 . 

16. A substance according to claim 15 and wherein said neutralizing comprises 



complementanly binding said KNA. 
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17. A substance according to claim 1 5 and wherein said neutralizing comprises 
immunologically neutralizing. 

18. A method for anti-viral treatment comprising neutralizing said RNA of claim 1 . 

19. Amethod according to claim 18 and wherein said neutralizing comprises; 

synthesizing a complementary nucleic acid molecule, a nucleic sequence of which 
complementary nucleic acid molecule is a partial inversed-reversed sequence of said RNA; and 

transfecting host cells witii said complementary nucleic acid molecule, thereby 
complementarily binding said RNA. 

20. A method according to claim 1 8 and wherein said neutralizing comprises 
iimmmologically neutraliziiig. 



51 



wo 2004/048S11 



PCT/IL2003/000998 




1/21 



wo 2004/048511 



PCT/IL2003/0OOS98 



FIG. 2 



riio 

BIOINFORMATIC 
GEISE DETECTION 
ENGINE 
TRAINING & 
VALIDATION 



SEQUENCED 
DNADATA 



PROTEIN 
FUNCnONDATA 



BIOINFORMATIC GENE DETECTION 
ENGINE 



NON-CODING GENOMIC SEQUENCE 
DETECTOR 



HAIRPINf 
DETECTOR 



DICER-CUT LOCATION 
DETECTOR 



TARGET-GENE BINDING-SITE 
DETECTOR 



FUNCTION & UTILITY 
ANALYZER 



BIOINFORMAHCALLY DETECTED 
GROUP OF NOVEL GENES 



2/21 



wo 2004/048511 



PCT/IL2003/000998 



FIG. 3 



^122 

TRAIN ENGINE ON KNOWN GENES 



HAIRPIN DETECTOR 
TRAINING & VALIDATION 



/-126 

DICER-CUT LOCATION DETECTOR 
TRAINING & VALIDATION 



/-128 

TARGET-GENE BINDlNG-SnEDETECrOR 
TRAINING & VALIDATION 





r /-130 


DETECT SAMPLE NOVEL GENES 




r /-132 


VALIDATE EXPREE 


.SION& FUNCTION 


OF SAMPLE GENES BY WET-LAB 



3/21 



wo 2004/04SS11 



PCT/IL2003/000998 



FIG.4A 



f EXPRESSH)^ fi 
RNADATA j 1 



^10 4 
SEQUENCED^ 
DNADATA I 



NON-PROTEIN CODING 
GENOMCSEQUENCEDiEIBCIOR 
mAJNING 



FIG. 4B 



NON-PROTEIN CODING 
GENOMIC SEQUENCE 
DETECTOR 



NON-PROTEIN CODING 
GENOMIC SEQUENCES 



GET KMA SEQUENCES 



SELECT 
NON-PROTEIN CODING 



LOCATE ON DNA 
(BLAST) 



EXPAND TO 
TRANSCRIPTION START & 
END SIGNALS 



DEFINITE EXPRESSEDNON- 
PROIEDSICQDING 
GENOMIC SEQUENCES 



GET DNA SEQUENCES 






PARSE 
NON-PROTEIN CODING 






LOCATE TRANSCRIPTION 
START & END SIGNALS 






PROBABLE EXPRESSED 

NQN-PR0IEINOODING 
GENOMIC SEQUENCES 



4/21 



wo 2004/048S1I 



PCT/IL2003/000998 



EIG.SA 



HAIRPnSf DETECTOR TKAINING 
& VALIDATION 



/-13 6 



NON-PROTEESf CODING 
GENOMIC SEQUENCES 





r /-I 14 


HAIRPIN I 


)ETECTOR 



c 



HAIRPINS ON GENOMIC 
SEQUENCES 



FIG.5B 



CALCULATE FOLDING & FREE- 
ENERGY 

T 

DETERMINE HAIRPIN FOLDING 

i 

DETERMINE IF SIMILAR TO 
KNOWN HAIRPIN-GENES 

i 

GROUP HAIRPINS INTO HAIRPIN 
'FAMILIES' 

i 

SELECT BEST 
HAIRPIN 'FAMILIES' 

r 

SELECT BEST HAIRPIN FROM EACH 
HAIRPIN 'FAMILY' 



5/21 



wo 2004/048511 



PCTAIL2003/000998 



FIG. 6A 



^13 8 



£ll 

DICER-CUT LOCATION 
DETECTOR TRAINING & 
VALIDATION 




DICER-CUT SEQUENCES 
FROMHAIRHNS 



6/21 



wo 2004/048S11 



PCT/IL2003/0CW998 



FIG. 6B 



GET DICER-CUT LOCATION OF KNOWN 
MICRO-RNA GENES 



REPRESENT 3' END 
(LOCATION, SEQUENCE, 
'BULGES') 



REPRESENT 5' END 
(LOCATION, SEQUENCE, 
^BULGES') 



TRAIN 
5' END RECOGNITION 
CLASSIFIER 



TRAIN 
3' END RECOGNITION 
CLASSIFIER 



DETECT DICER-CUT LOCATION 
OF SAMPLE NOVEL GENES 



VAUDATB BY WET LAB 



7/21 



wo 2004/048511 



PCT/IL20037000998 



FIG. 6C 



GET 

NOVEL VGAM PRECURSOR 




mi 



wo 2004/048S11 



EIG.7A 



PCT/IL2003/000998 

-142 



128 



DICER-CUT SEQUENCES 
FROM HAIRPINS 



TARGET-GENE 
BINDING-SITE 
DETECTOR 
TRAINING & 
VALIDATION 



[CES \ 





POTENTIAL HOST 
TARGET GENE 



TARGET-GENE BINDING-STTE 
DETECTOR 




^144 



POTENTIAL NOVEL TARGET- 
GENES 
HAVING BINDING-SITE/S 



FIG.7B 



J 



BLAST COMPARE 



FILTER BLAST RESULTS 



EXPAND BINDING SITE 



CALCULATE FREE-ENERGY & 
SPATIAL STRUCTURE 



ANALYZE 
SPATIAL STRUCTURE 



wo 2004/048S11 



PCT/1L2003/000998 




10/21 



wo 2004/048511 



PCT/IL2003/000998 



FIG. 9 



VGR PRECURSOR KNA 



II If It YGRFOIDED 

ffiip! 

I \ 4 
w% it i% 



VGAM3 

"rna 



^ TARGET RNA 
SITE _ 

r Tojyve-TARGET" ! 
PROTEIN 



- BINDING ^ , V- VGAM2- 
SITE ^ TARGET RNA 



VGAM2-TARGET 
PROTEIN 



T 



VGAMl 

RNA 



BINDING "^i!^^ VGAMl ■ 
TARGEl 



TARGET RNA 



VGAMl-TARGET ' 
PROTEIN 



11/21 



wo 2004/048511 



PCT/IL2003/000998 



MG. 10 



unLmES 



VGAMGENE 
DETECTION 



ANn-VGAM 
GENE 
THERAPY 



12/21 



wo 2004/048511 



PCT/IL2003/000998 



FIG. IIA 

■VGAMGENE 



\ 



VGAM PRECURSOR 
/HRNA 



-DICER COMPLEX 



VGAM FOLDED 
,rPRECURSORRNA 

T^VGAMBNA 



-VGAMRNA VGAM 

/^TARGET RNA 



V BINDING i 



VGAM TARGET 
PROTEIN 



FIG. IIB 



I 



VGAM PRECURSOR 
/HRNA 



DICER COMPLEX 



VGAM FOLDED 
"PRECURSOR RNA 



-VGAMRNA 

-ANItVGAMRNA 

VGAM 
-^TARGET 
RNA 



[binding 




saiE 






VGAM TARGET 






PROTEIN 





13/21 



wo 2004/048511 



PCT/IL2003/000998 



5 « § 

o < p 

z o i 



s. in 
o S 
c < 

"11 
55 



1 5 2 

0- o 



■2 I 

I i 1 



1 

5 



SUBSTITUTE SHEET (RULE 26) 



wo 2004/048511 



PCT/IL2003/000998 



1 2 3 4 5 



FIG. 13 



Library - + -+-+- + - + 

6 7 8 9 10 11 12 13 




Library 



Library - + - + - + - 

10 11 12 13 



Library . + _ + - + 

c a 

Library - + 



D 

Library - + 



15/21 



SUBSTITUTE SHEET {RULE 26) 



wo 2004/048511 



PCT/IL2003/000998 




SUBSTITUTE SHEET (RULE 26) 



wo 2004/048511 

EST72223f705 nU 



PCT/IL2003/000998 



FIG. 15A 



EST72223 sequence: 

rf^nTTATTAnAnf^ATTrTfinTnATnnnARRG TGAGGTAGTAAGTTGTATTG 

HGTGGGGTAGGGATATTAGGCCCCAATTAGAAGATAACTATACAACT MIR98 

TACTACTTTCCCTGGTGTGTGGCATATTCACACTTAGTCTTAGCAGTGTTGCC 

TCCATCAGACAAAGTTGTAGATGTTCCTTGGATAATTTGGACTGGAAGAAAAGA 

GACATGGAAGGGGACAGATGGTGTTTAGGGT GA6G CAGATGTCATTATAAAGT 

GACTTGTCTTTCATTAATTGGAGCATATAATTA TTTT ACCTTTGGGCATGAACTC 

ATTTTGCTATTCTTCAACTGTGTAATGATTGCATTT TATTA GTAATAGAACAGGA 

ATGTGTGCAAGGGAATGGAAAG CATACT TTAAGAATTTTGGGCCAGGC6CGGT 

GGTTCATGCCTGTAATCCCAGCATTTTTGGGAGGCCGAGGCGGGTGGATCAC 

CTGAGGTCAGGAGTTCGAGACCAACCTGGCCAACACGGCGAAACCCCGCCTC 

TACTCAAATACAAAAATTAGCCAGGCTTGGTGACACTCGCCT6TGGTCCCAGC „ 

TArTnAPiriARrinTGAGGCAGGA QAATTGCTTGAACCCAGGAAGTGGAG GAM25 

GCTTCAGTGAGCTGAGAACACGCCACTGCACTCCAGTCCTGGGCAAC 

AGAGCAAGACTCTGTCTCAGGAAAAAAAAAG 
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5'UTR SEQUENCE (5* TO 3') OF HIV-1(U5-R) 

G GTCTCTCTGGTTAGACCAGATCT GAGCCTGGGAGCTCTCTGGCTAACT 
AGGGAACCCACTGCTTAAGCCTCAATAAAGCTTGCCTTGAGTGCTTCAAGTA 
GTGTGTGCCCGTCTGTTGTGTGACTCTGGTAACTAGAGATCCCTCAGACCCTT 
TTAGTCAGTGTGGAAAATCTCTAGCAGTGGCGCCCGAACAGGGACCTGAAAG 
CGAAAGGGAAACCAGAGGAGCTCTCTCGACGCAGGACTCGGCTTGCTGAA 
GrGrGTArG GCAAGAGGrGAGGGGCGG CGACTGGTGAGTACGCCAAAAA 
TTTTGACTAGCGGAGGCTAGAAGGAGAGAG 
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